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AN APPROXIMATE SOLUTION OF THE EQUATIONS OF MOTION 
FOR ARBITRARY ROTATING SPACECRAFT* 

By Peter R. Kurzhals 
Langley Research Center 

SUMMARY 

The assumption of small changes in the inertia parameters has been used to derive 
approximate rotational equations of motion for arbitrary spinning spacecraft in the small 
angle and rate regime. Complex representations are introduced to define the rate and 
attitude errors produced by applied disturbances, and analytic solutions are obtained for 
the steady spinning mode and for the spin-up and despin mode. 

Solutions for the steady spinning mode consider both the uncontrolled and the con- 
trolled spacecraft motion for characteristic disturbances. These disturbances include 
initial errors, externally applied torques, and instantaneous and periodic mass motions 
within the spacecraft. The errors induced by the disturbances are described by the 
error component time histories and by vector traces of the complex error representa- 
tions. Upper bounds of the errors are developed for the uncontrolled motion, and the 
required control techniques and control systems are examined for the controlled motion. 

Solutions for the spin-up and despin mode consider extensible spacecraft modules 
connected by struts or cables. Fuel consumption relations are derived for several 
extension techniques, and optimization of the extension techniques is shown to yield 
appreciable fuel savings. 

Comparisons of the analytical solutions and exact solutions obtained by numerical 
integration of the complete equations of motion are found to yield excellent agreement, 
and the applications of the approximate solution are illustrated for a manned orbital 
research laboratory and a large spinning space station. 


*The information presented herein was submitted as a thesis in partial fulfillment 
of the requirements for the degree of Doctor of Philosophy in Aerospace Engineering, 
Virginia Polytechnic Institute, Blacksburg, Virginia, May 1966. 


INTRODUCTION 


Proposed spacecraft, such as the manned orbital laboratory (ref. 1) and manned 
interplanetary vehicles (ref. 2), may use rotation about a maximum axis of inertia to 
provide spin stabilization and to produce artificial gravity for the crew. These space- 
craft will be subjected to variable torques arising from both internal and external 
sources (ref. 3) and will undergo wobbling motions as a result of these torques. Since 
the wobbling motions (ref. 4) produce attitude errors (which may affect the spacecraft's 
power system and experiments) and oscillatory rates (which may lead to discomfort and 
nausea of the crew), an analysis is required to determine the magnitude of any such atti- 
tude errors and body rates for the spacecraft under consideration. 

In order to carry out this analysis, the spacecraft's equations of motion with 
varying inertias and torques must be integrated to define the spacecraft response for the 
anticipated applied disturbances. In the past, such a solution has required high-speed 
computing equipment for the numerical integration of the equations of motion and has 
consumed a large amount of computer time to assess the effects of a range of disturb- 
ances for a particular vehicle configuration. 

Because of the rather limited application of these results, an approximate analytical 
solution of the spacecraft's equations of motion would be of considerable value. The 
closed-form solution could be used to determine attitude errors and body rates intro- 
duced by "worst case" type of disturbances and would define instability trends that might 
result from applied torques. In addition, such a solution would allow a direct evaluation 
of the effects of changes in both the spacecraft's configuration and the disturbances on the 
spacecraft's motion. 

Approximate analytical solutions of the equations of motion for an arbitrary rotating 
spacecraft may be obtained for linearized governing equations and Euler angle transfor- 
mations. Several such solutions have been obtained for the simplified equations of motion 
corresponding to symmetric or near- symmetric spacecraft. Leon (ref. 5) and Thomson 
(ref. 6) have developed attitude and rate relations for spinning near- symmetric bodies by 
considering a vectorial representation of the total errors. Thomson and Fung (ref. 7) 
have also investigated the stability of near- symmetric spinning space stations and have 
defined regions of instability for an example vehicle. In addition, Poli (ref. 8), Buglia 
and associates (ref. 9), and Loebel (ref. 10) have derived expressions for the attitude and 
rate histories of symmetric spacecraft by linearizing the equations of motion. 

Several analytical solutions for a nonsymmetric spinning body with constant iner- 
tias have also been obtained. Exact solutions for a torque-free body were developed by 
Routh (ref. 11) and MacMillan (ref. 12) in terms of Poinsot’s construction and elliptic 
functions and by Whitbeck (ref. 13) in terms of a phase plane approach. An approximate 
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method which shows good agreement between the nonlinear and linearized results for a 
vehicle under applied torques was presented by Suddath (ref. 14). 

Various other analytical approximations are discussed in the literature (refs. 15 
to 20). Existing solutions, however, have considered either very special cases of non- 
symmetric spacecraft or have been restricted to particular symmetric or near- 
symmetric spacecraft with specified disturbances. Such results cannot be applied to the 
general case of a nonsymmetric spacecraft with varying products of inertia and applied 
torques and offer little information on the properties of the motion of such spacecraft. 
Furthermore, the form of these solutions has made the determination of upper limits for 
the total attitude and rate errors difficult since the amount of computational time required 
to define the error boundaries is, in general, prohibitive. 

The present analysis develops a solution technique for arbitrary rotating spacecraft 
with variable disturbance functions. The complete equations of motion for nonsymmetric 
vehicles are linearized and solved with time-varying forcing functions and products of 
inertia. General and particular solution functions are determined and are used to gen- 
erate rate and attitude expressions corresponding to the variable forcing functions. A 
complex vector representation is introduced to define both error time histories in com- 
ponent form and the total angular and rate errors. 

A number of disturbances are considered for both nonsymmetric and symmetric 
spacecraft, and the corresponding solutions are examined for the uncontrolled and con- 
trolled cases. Upper bounds of the total errors are defined and body-fixed and inertial 
traces of the total errors are analyzed. A method of selecting control commands is also 
presented. 

SYMBOLS 

Aj complex coefficients of forcing function (see eqs. (Bl) and (30)) 

a power conversion factor (eq. (254)) 

a,b,c,d,e characteristic coordinates for error traces 
Bj complex coefficients of total rate error relation (eq. (43)) 

Cj complex coefficients of total angular error relation (eq. (44)) 

d 1> d 2’ D 3’ D 4’ D 5 spin-up parameters (see eqs. (295) and (297)) 


3 



Eri.Ei.ErwEo.E^ forcing function coefficients for variable inertia products (see 

table IV and eq. (112)) 

EgjEgjEiy solution function coefficients for variable inertia products (see eqs. (113) 
and (115)) 

F complex forcing function (eq. (30)) 

F complex solution function for body rate errors 

F complex solution function for Euler angle errors 

forcing functions for the controlled spacecraft (see eqs. (213) and (214)) 

G complex actuator torque 

Gg actuator stall torque 

g complex control torque (eq. (209)) 

H rigid-body angular momentum 

I moment or product of inertia 

IpQ spacecraft moment of inertia without moving particles 

I rz radial inertia product, I xz + iIy Z 

Igp specific impulse 

I T total impulse (eq. (284)) 

Kj x ,Kjy control gain for X- and Y-axis, respectively (eq. (209)) 

kj x ,kjy nondimensional control gain for X- and Y-axis, respectively (eqs. (210)) 

L nondimensional external moment (see eqs. (A30)) 

l distance between mass centers of manned and counterweight modules 
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effective spin-up moment arm (see eqs. (294)) 


external moment, M x + iMy 

spin-up moment 

mass of moving particle 

mass of spacecraft and moving particles 

order of [a] 

external force or thrust 

frequency for moving particle 

constant defining center-of-mass change (eq. (A25)) 

position vector of moving particle measured with respect to origin of fixed 
coordinate system 

position vector of origin of spacecraft coordinate system measured relative 
to origin of fixed coordinate system 

damping ratio 

position vector of moving particle measured with respect to spacecraft 
coordinate system 

position vector of composite mass center measured with respect to space- 
craft coordinate system 

Laplace transform variable 

signum function, denoting sign of characteristic coordinate or variable 

constant step or impulsive torque, T x + iT y 

cross coupling torque applied to spacecraft by control system 


t 

U(t) 

V 

w 

X,Y,Z 

x,y,z 

u,v,w 


Q! 


"I 

* 


y i 

6(t) 

A 


€ 






time 

time constant 

unit step function (eqs. (B4)) 
effective forcing function 
weight 

reference axes 

scalar component of r" along X-, Y-, and Z-axis, respectively 

nondimensional scalar component of r" along X-, Y-, and Z-axis, 
respectively (see eqs. (A30)) 

complex angular position error, cp + id (see fig. 4) 

limiting gyro gimbal angle (eq. (272)) 

complex inertial position error (eq. (14)) 

argument of Bj 

argument of Cj 

impulse function (eqs. (B5)) 

quantity denoting magnitude order of cl)p, 8cp, 6^, Lp, Pp, and non- 

dimensional inertia terms for moving particles 

nondimensional moment or product of inertia (see eqs. (A29)) 

modified Euler angles (fig. 2) 

precession rate parameters, positive when I z is maximum inertia and 
negative when I z is minimum inertia (eqs. (24)) 
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nondimensional mass (see eqs. (A29)) 

nondimensional force (see eqs. (A30)) 

constant positive spin rate (eq. (20)) 

nondimensional time (see eqs. (A29)) 

angular coordinate used in total error traces 

total angular rate vector of spacecraft axis system 

scalar component of S2 along X-, Y-, and Z-axis, respectively 

complex rate error, S2 X + if2y (fig. 3) 

damped natural frequency 

synchronous wheel speed 

nondimensional scalar component of £2 along X-, Y-, and Z-axis, 
respectively 


Subscripts: 


a,b,c,d value for corresponding characteristic coordinate 
c counterweight module 

CM extension at constant momentum 

CR extension with constant spin rate 

CT extension with continuous thrust 

d disturbance 

e value after extension 


F 


fixed coordinates 



final value after spin-up 


control moment gyro 
gyro gimbal 

intermediate coordinates 

intermediate value before extension 

value for jth term or mass where j = 1, 2, 

component for X- or Y-axis with k =£ l 

upper bound 

momentum 

manned module 

maximum value 

summed value for moving particles 
value referred to origin of spacecraft body 
initial value 
power 

component for X-, Y-, or Z-axis 
component for XY-, XZ-, or YZ-plane 
reaction control 


residual value 



SU spin- up fuel 

s value for spacecraft mass center 

T total 

x,y,z component for X-, Y-, and Z-axis, respectively 

xy, xz, yz component for XY-, XZ-, and YZ-plane, respectively 
W reaction wheel 

Mathematical notation: 

. over a symbol denotes derivative with respect to time 

— over a symbol denotes a vector 

~ over a symbol denotes general solution function corresponding to initial 

rate and attitude errors 

with a symbol denotes particular elements of a vector component along 
spacecraft axis 

| | absolute value 

The quadrant for the angles corresponding to the inverse trigonometric functions 
tan” ^["3 is determined by the sign of the numerator and denominator of the term within 
the brackets. When both numerator and denominator are positive, the angle is in the 
first quadrant; when the numerator is positive and the denominator is negative, the angle 
falls in the second quadrant; when both numerator and denominator are negative, the 
angle is in the third quadrant; and when the numerator is negative and the denominator is 
positive, the angle falls in the fourth quadrant. 

All square root terms in this analysis are principal, positive values. These values 
may be positive real or positive imaginary numbers. 
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The use of logarithmic figures in this report makes the incorporation of both U.S. 
Customary Units and the International System of Units (SI) impractical. U.S. Customary 
Units are accordingly used in the report text, figures, and tables. Conversion factors to 
SI Units are as follows: 

Quantity 

H 

W 

I 

X T 

T,G 

Q 

m 

I 

x,y,z 

a 

Kl 

k 2 ,k 3 

d 2 

D3 

i 

PROBLEM FORMULATION 


U.S. Customary Unit SI Unit 


ft-lbf-sec 

1.356 N-m-sec 

lbm 

0.454 kg 

slug- ft 2 

1.356 kg-m 2 

lbf-sec 

4.45 N-sec 

ft-lbf 

1.356 N-m 

slug 

14.59 kg 

slug 

14.59 kg 

ft 

0.3048 m 

ft 

0.3048 m 

lbm/watt 

0.454 kg/watt 

ft-lbf-sec/ rad 

1.356 N-m-sec/rad 

ft-lbf/ rad 

1.356 N-m/rad 

lbf 

4.45 N 

ft 

0.3048 m 

ft/ sec 

0.3048 m/ sec 


Spacecraft Motion 

The rotating spacecraft will be related to the reference system shown in figure 1. 

A set of X,Y, Z axes fixed to the spacecraft is used to describe the rotational motion 
of the spacecraft with respect to a set of intermediate Xj,Yj, Zj axes. The intermediate 
axes translate without rotation in inertial space but always remain parallel to a set of 
Xp,Yp,Z F axes fixed in inertial space. 

The inertial attitude of the spacecraft may be defined by means of three modified 
Euler angles which determine the relative motion between the X,Y, Z and Xj,Yj,Zj 
axis systems. These modified Euler angles, as illustrated in figure 2, result from three 
consecutive rotations. The first rotation, about the Zj-axis, carries the Xj- and Yj-axes 
through an angle xp measured in a horizontal plane. The second rotation, about the 
new Yr axis, then takes the Xj- and Zj-axes through an angle 6 measured in a vertical 
plane. Finally, the third rotation, about the new Xj-axis, carries the Yj- and Zj-axes 
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Figure 1.- Reference system for rotating spacecraft. 



Figure 2.- Vector transformation between spacecraft axes 
and intermediate reference system. 


through the angle cp measured in an 
inclined plane to give the X-, Y-, and Z-axes. 

The modified Euler angles can be 
determined by expressing the rotations if/ 

9, and ip in terms of the angular rates O^ 
Oy, and O z about the vehicle axes. The 
vehicle angular rates then can be found from 
a solution of the vehicle moment and force 
equations. The resultant expressions for 
Ox, Oy, and O z are substituted into the 
Euler angle transformations, which now 
reduce to differential equations in if/, 9, 

(p, and t. The solutions of these equations 
give the attitude of the spacecraft relative 
to the intermediate axes and thus determine 
the angular motion of the spacecraft. 

Assumptions 

To make the general nonlinear equa- 
tions of motion amenable to analytical treat- 
ment, the assumptions 

Ox O z 

Oy « O z 

T > (1) 
sin 9 = tan 0=6 cos 8 = If 

sin cp = tan cp = cp cos cp = 1^ 

were introduced in the moment, force, and 
Euler angle relations developed in appen- 
dix A. The further assumption that the non- 
dimensional inertia terms associated with 
any mass particles moving with respect to 
the spacecraft were small was also made to 
linearize the equations. The resultant 
method of reduction to linear form and the 
range and validity of these assumptions are 
discussed in appendix A. 
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Governing Equations 

With the assumptions of the preceding section, the equations of motion reduce to 



_ 


f 

n 

fix + 

(rrH 

Hi -? 5 

II 

dr 

M x + ^z(ixz “ lyz^z) + 

Z m j( z 3 x j ■ x j z i) “ m s( z s x s - x s z s ) 
3=1 


n 


+ 1 - y fj) ■ m s( z s y s - y s z s)J 

j=i 


S2 n 


( 2 ) 


% 



r 

[W-] 

*> 

II 

-ri- 

< 


Y mjfZjYj - yjzj) - m s (z s y s - y a z s ) 


n 


+ I m i( x j z i " z 3 x i) " m sK z 
J=i 


S - z s x s 




fi„ 


(3) 


^ * T z ^ + 1 m i( y j*j ■ x s y j) - m s(y s *s - x sy s ) 

J=i 


(4) 


The inertia terms are 


n 


2 , „ 2 ' 


J x = X xO + X m i( y jo 2 + z io 2 ) - m s^y so “ + z so‘ 

3=1 

n 

h = 1 yO + J J m j( x jo 2 + z jo 2 ) - m s( x so 2 + z so 2 

3=1 

n 

iz = ! z o + £ m j( x j 2 + yj 2 ) - m s( x s 2 + ys 2 '' ^ 

3=1 


n 

*xz = ^ m j x j z j ‘ m s x 


s z s 


3=1 

n 


i yz = ^ m jyj z i - m s y s z s 

3=1 


j 


(5) 
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and, consistent with the assumptions, the moments of inertia are taken to be constant in 
equations (2) and (3) but are allowed to vary in equation (4). The associated inertia 
derivatives then become 


ly = 2 


n 


-T\ 


X m j( x i k 3 + y fj) “ m s( x s x s + ysYs) 

i=l 

n 

ixz = ^ m i( x j^3 + z i x j) “ m s( x s z s + z s x s) ^ 


1=1 

n 


iyz = J m j(yj z i + z jyj) - m s(y s z s + z s y s ) 

i=i 


j 


( 6 ) 


where xj,y-j,Zj denote the position coordinates of the mass mj moving with respect to 
the spacecraft, and 


x 


s 


y — 

Zy m S 3 


y s = 


n 



m i 

y 3 


(7) 


z s 


n 

Y^ Zi 

l_j m s ] 


denote the position coordinates of the composite mass center for the spacecraft and the 
moving masses. 

The spin rate S2 Z is obtained by integrating equation (4) as follows: 


°z = pko^zo + J m z dt + J £ m j(yj x j - x fj) - m s(y s x s - x sy s ) 

z l k=l 


( 8 ) 


where the first term in the braces represents the system's initial angular momentum and 
the remaining terms account for changes in the spin rate due to the applied torque M z 
and to the accelerations of the moving masses. 

Solutions to the spacecraft equations of motion may be obtained by first deter- 
mining S2 Z from equation (8) and then integrating equations (2) and (3) simultaneously 
to find Six an d Oy. 


13 



These solutions can be substituted into the linearized Euler angle transformations 


cp — £2x + fi^ 0 

(9) 

0 = fiy — ^2 Cp 

(10) 

II 

S3 

N 

(11) 


so that 


i r ~ n 

* = J ^^zo^zo + J m z dt + J - x jyj) - m s(y s x s - x s y s ) 

l ^ 


dt)dt (12) 


The solution of equations (2), (3), (9), and (10) then defines the motion of the spacecraft 
in terms of the time histories of the body rates and Euler angles. 


Total Errors 


The body rates fix and fiy are the undesired rate components produced by the 
applied disturbances and are referred to as the rate error components. Similarly, the 
Euler angles <p and 6 describe the unwanted attitude deviations that result from the 
application of the disturbances and are referred to as the attitude error components. 

The solutions for both rate and attitude error components follow directly from the pre- 
ceding section and are found as time dependent components along the body and inertial 
axes. 

In practice, the total errors are of primary concern. For example, the time vari- 
ation and maximum value of the total angular velocity error in body-fixed coordinates 
must be known to assess possible crew discomfort due to wobbling motions. The time 
variation and maximum value of the total angular position error with respect to inertial 
space is needed to determine possible effects on the spacecraft experiments and power 
system. The effects of removal of a disturbance on the residual spacecraft motion are 
also of interest. 


Both the total angular position and the total body rate errors may be developed by 
using a complex vector representation (ref. 5). The total angular rate error fixy can 
be obtained by vector addition of the body rates fix and fiy, as shown in figure 3. 
Mathematically, fixy may be written as 


fixy = fix + ifiy = 





(13) 


Similarly, the total attitude error a in body-fixed coordinates can be considered to be 
the vector sum of the small Euler angles cp and 6, as illustrated in figure 4. To 
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Y 


Figure 3.- Vectorial representation of total angular rate error. Figure 4.- Pseudovectorial representation of total attitude 

error. 



transform this pseudovector to the intermediate coordinate system, the body coordinate 
system must be rotated through the angle The total inertial error a j is then 


= ae^ = ((?+ id)e^ = \Jcp^ + 6^e 


,rp + tan^d) 


Physically, represents the trace of the Z-axis projected on the XjYj-plane and 
f2 xy represents the trace of the total rate error vector in the body-fixed XY-plane. 

Since differentiation of equation (14) yields 

At — (At (1 fSi 


and noting that 


otj. = (a + i^a)e 


of + — Oxy 


from equations (9) and (10), equation (11) may be used to develop the relation 


. / • \ iO^t iilgt 

= [a + iOt^aje = Slxye 


The magnitude of the rate of change of the inertial attitude error is thus equal to the 
magnitude of the rate error for the small angle regime. 

Integration of equation (17) results in 


Oj = + \ Six ye ^ dt 
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as the solution for the inertial attitude error vector. The attitude error in body coor- 
dinates becomes 

a = a-je' 1 ^ = a^e"^ 2 * + e"^ 2 * J Si^e^dt (19) 

and both the attitude errors oij and a can be directly developed from the rate error 
expression. 

If only total error vectors are desired, the time solution of the equations of motion 
for f2xy may be followed by application of equations (18) and (19) to yield a j and a. 

If the rate and Euler angle components are of interest, the direct solution of the linear 
differential equations (2), (3), (9), and (10) is preferable. 

In the present analysis, solutions were first developed in the form of time histories 
for the error components. The transition to the total error form was then made by sub- 
stitution of the vector expression for the resultant rate error in equations (18) and (19). 

Solution Approach 

The solutions of the equations of motion for arbitrary rotating spacecraft can, in 
general, be divided into two types: namely, those associated with the spin-up and despin 
modes and those associated with the steady spinning mode. 

The spin-up and despin modes may involve the extension and retraction of cable- 
connected counterweight modules and thus could produce major and rapid changes in the 
moments of inertia for the spacecraft. During these modes other disturbances, such as 
crew motions and applied torques, will necessarily be restricted and only the solution for 
the spin rate and angle, as given by equations (8) and (12), need be considered. The effi- 
ciency of various spin-up and despin methods using constant spin rate, constant cable ten- 
sion, or similar schemes can be readily evaluated from these equations. 

For the steady spinning mode, the variations of all total moments of inertia due to 
the moving masses associated with a particular crew motion are small in comparison 
with the constant spacecraft inertias I^q, IyQ, and I z q. The assumption that the total 
moments of inertia and I z retain their initial values throughout the crew 

motion (see appendix A) is made for this mode. 

In addition, disturbance moments due to crew motions and applied torques now act 
primarily about the spacecraft X- and Y-axes. Any torques about the spacecraft Z-axis 
can be neglected during a particular disturbance since the resultant change (refs. 8, 14, 
and 19) in the spin rate will be small in comparison with the initial spin rate. 

In accordance with these assumptions, the spin rate may be approximated by its 
constant value at the initiation of a particular disturbance 
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( 20 ) 


£2 Z — 


*zQ^zo ; a 

I z 


for the evaluation of the effects of that disturbance on the spacecraft motion in the steady 
spinning mode. The value a can be taken to be positive without loss of generality. 

Since the spinning mode occurs for the major portion of the spacecraft lifetime, 
this mode is first analyzed in detail and several spin-up and despin techniques are then 
considered. 


ANALYSIS OF SPINNING MODE 

For the spinning mode, the moments of inertia take on their initial values imme- 
diately after initiation of the disturbance and remain constant for the duration of the dis- 
turbance. The inertias may thus be computed from equations (5) as 

n '' 

I x = l xo + J m j(yjo 2 + z jo 2 ) - m s(y so 2 + z so 2 ) 


; jQ 2 ) “ m s( x S0 2 + Z SO ] 


j=l 

n 

ly = IyO + ^ m j( x jo“ + z jO 

j=l 

l z = x zo + £ m j( x jo 2 + yjo 2 ) - m s( x so 2 + y so 2 ) 
i=i 




( 21 ) 


J 


and the governing equations may now be developed directly from equations (2) and (3). 

Equations of Motion 
The equations of motion reduce to 


S2x + 


m j( Z A “ x j*^j) " m s( z s* x s “ x s z s) 

j=l 


( % = £ j** ' My - [ iyz<7 + ixzo2 ] + E xza - iyza2 ] 

n 

+ £ m j( z j*y j + zjVj - - yfi) - m s( z s*y s + z sv s - y s * s - y s z s) 

5=1 

- [(^V 5L ) a m i( z ji r i - y 5 fe j) - m s( z s y s - y s z s)| <j + f m i( x j z j - z i x j) - m s( x s z s - z s x s) 


( 22 ) 
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and 




% " ) M y + 


(^5~)° r M x + (^5“) CT (^ ~ V* 0 ^) + ^"yz° r + 

+ 2 m j( x j z j + X j Z i " Z j' X j ■ Z j X j) ■ m s( X s’ Z S + X S Z S - z s’ X S - z#s) 


Z m i( z i y 3 ' y f j) ■ m s(vs - y s z s) 
13=1 


fe) ; 


'j( z 3 x 3 - x j z j) - m s( z s x s - x s z s)>°+ X m j( z j y j • y j z i) - m s( z sy s - ys'zs) 

li=i J 3=1 


( 23 ) 


For simplicity of notation, introduce the precession rate parameters 




^Z “ lx 


lx 






X^ — X X Xy 


(24) 


so that equations (22) and (23) become 

n 

fez - iyz( a+ V) - I xz°^y| + ^ m j f( a + X y)( z f j " x fj) 

j=l L 

» y ( Z) y) - yjij) + (zj'y j + Vi - yj'z'j - yjSj)J - m s |(a + x y )(z s x s - x s z s ) 


+ X = ” ^y M y 


- aXy^z s y s - Ys^s) + ( z s Y s + ^sYs " Ys z s Ys z s) 


and 


£2y + X f2y — + X x M x + a 1 


n 

; + ixz( a + ^x) “ lyz^xl + ^ m j[( CT + ^x)(zjyj “ YjZj) 

i=i 

+ aXx^ZjXj - XjZj) + (Xj-Zj + XjZj - z-x j - z-x^ - m s |(ff + A x )(z s y s - y s z s ) 


+ aX X ^Z g Xg - XgZg^ + (x S Z g + XgZg Z g X g ZgXg^ 
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or 


and 


where 


£^ x “J* A — F x 


^2y + X £2y — Fy 


(25) 

(26) 


F x = f<M x -X y My + 


p- — » n 

T’xz - iyz( a + A y) - Ixz^y) + £ m > f( CT + X y)( z i x j * x i z i) “ CTX y( z jyj - yj z j) 

3=1 

+ (zj-yj + Zjyj - yj'z'j - yjZjjj - m s |(o- + A y )(z s X s - x s i^ - cxA y (z s y s - y s z s ) + (z s 'y s + z s y s - y s 'z s - y s ‘i s )jj> 

= g^x - x y M y + I “jfjfy'i + ( 2ff + X y) x j - a {° + 2 x y)yj - + z j{^j + 2oi j - - *j{*j + Vij 

■ z jyj] - m s[ z s{y s + ( 2or + v)^ s - + 2x y)ys - °^^y x sj + z sj^s + 2°x s - ^yj} - z s j^s + V x ®} ~ z s^, 


(27) 


and 


y - i^y + X x M x + a [^yz + 3xz( CT + *x) ' Iyz°*x) + ^ m jJ{ CT + X x)( z jyj - yj z j) + aX x( z jXj - XjZj) 

+ (xjZ*j + XjZj - Zj-Xj - ZjXj)] - m s |ja + A x )(z s y s - y s z^) + aA x (z s x s - x s z s ) + (x s ‘z s + x s z s - z s x‘ s - z s x s || 

J m i[ z if*j - ( 2ct + x *)yj - 2x x) x i + ^xyj) + z j{^j - 2 °yj - Aj) - Sj{*j - x x yj} - -z^jj 

[zs^x's - (2<J+ X x )y s - a(cr+ 2X x )5% + a 2 X x y^ + z s js s - 2oy s - Agj - z s |k s - X x yJJ - 'z s x s | 


” ^\ M y + x x M x 


Adding equations (25) and (26) in quadrature and referring to equation (13) yields 


with 


The solution of equation (29) is 


£X^y -|- X £2 X y — F 


F = F x + iFy 


^^xvo — 

= — — ■ sin Xt + f2xyo ^os Xt + F 


(28) 

(29) 

(30) 

(31) 


where F is obtained by replacing functions of t in F by the corresponding particular 
solution functions given in appendix B. 
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The initial conditions at t = 0 are 


f^xyo - ^xo + i^yo 


A y I y 


S2xy0 _ " \TTr y° + 1 Ur 0x0 


A x^x 




( 32 ) 


The particular contributions of an applied disturbance to the initial errors are included 
in the Laplace formulation of the solution terms. 

Substitution of equations (32) into equation (31) then results in expressions for the 
total error flxy and its components fix and Oy» The spin rate S2 Z is found from 


Sly — 


^zO^zo 


= a 


and all body rates have thus been defined. 

The Euler angle differential equations can be written as 

cp + cfi(p = fix + CTflyl 

0 + a^0 = fly - oflxj 

and after adding in quadrature 

of + o 2 of ~ flxy —iofl^y 

Substitution for the rate error in equation (34) yields 


.. o 

Of + CT^O! = 


which has the solution 


a, 


iafl. 


Xflxyo J 


xyo 


sin Xt + 


^xyo _ i a ^xyo 


cos Xt + F - iaF 


a - — sin at + ofQ cos at + — *•< 

u - X" 5 


- - 

i) 

^xyo " * a ^xyo 

cos Xt - cos at ' 


*%yo - ^ 


+ F - iaF 


sin Xt - - sin at 


(33) 


(34) 


(35) 


(36) 


where F and F are obtained by replacing functions of t in F by the corresponding 
particular solution functions of appendix B. 
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Initial conditions at t = 0 are 


= <Pq + i0 o 

do = ^xyo - i GQ b[ 


and the Euler angle i// can be determined from the relation 

i// = at 

This relation completes the development of the Euler angles. 


(37) 


(38) 


Since the terms involving the initial errors have the same form for all disturbances, 
introduce 

~ ~ /v o <5 X VO 

+ i£2y = sin At + S^xyo cos At 


— £2 X iOy 


s * <f> + - 


A r^xyo _ ia£2xyo 


a 2 - X 2 


>sin at + ( Oq 


^xyo ~ ^ a ^xyo 


a^ - A^ 


cos at\ (39) 


, i f r 

o2 - X 2 \L 

or in component form 


iafi- 


X^xyo j- 


'xyo 


sin At + 


A V I 


f^xyo ” ia^xyo At 


Q x = ^2x 0 cos At - [-XZjftyo sin At 




Sly = Sly Q COS At + [-^r-'J^xo sin ^ 


and 


- [*> - ^ 
= [* 


cos at + 


Ix^XO 


6r\ + T 

° al 7 


sin at + 


(40) 


I y%o 


al 7 


cos Xt + 


■^x^x^xol 


aXI, 




sin Xt 


Ix^xo 

al z _ 


cos at 


<Po ~ 


IyS2y 0 


al, 


z J 


sin at - 


*x^xo 
L ffI zJ 


cos Xt + 


V x y^yo 

aAI 7 


} (41) 


sin Xt 


All coefficients in these equations may be evaluated from the initial conditions £2xo, 
£2yo, a, (p Q , e 0 , and the initial moments of inertia. 
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The terms involving the applied disturbances can be similarly put into component 
form, so that the body rates and Euler angles can be found by equating real and imaginary 
parts in 

f2xy = fix = flxy F 

O!=<p + i0 = 2+ F- icrF 

where F, F, and F are taken from appendix B. 

Total Errors 

The total rate and attitude errors may be put into a somewhat simpler form by 
expressing all trigonometric terms occurring in these errors in exponential form. Thus, 

^ = I I (43) 

j=l h=0 

where j and h vary over a finite range of integers. The complex constants Bj and 
the real constants must be evaluated for a particular disturbance. 

The total angular error in inertial space, as defined by equation (14), can be sim- 
ilarly expressed as 

u v t 

= (<p + i 0)e i ^ = ^ ^ Cjt h e ^ (44) 

j=l h=0 

where j and h again remain finite integers with Cj and yj determined for a 
specified disturbance. 

In engineering applications, the maximum magnitudes of these errors are also of 
interest. Since the exact solution for the maximum error magnitudes requires an itera- 
tive determination of the zeros of the magnitude derivative — and this requires consider- 
able computing time — an alternate method of defining upper bounds for the errors is 
preferable from the practical standpoint. It is noted from equation (43) that 
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|°xy| ~ 


x ! -a 1 * 

j=l h=0 


f g 




j=l h=0 

and similarly from equation (44) that 


*XXh 

j=l h=0 

s X X NM * Klu 


lim 


( 45 ) 


ii 

jr 

K 

jr 

X X 


j=l h=0 


- X X l c «* 

3=1 h=0 


u 


s 2 X Nl th 

i=i h=o 


= O' 


lim 


(46) 


These upper bounds provide limiting values of the error magnitudes which are adequate 
for assessing the effects of particular disturbances. More accurate estimates of the 
absolute maximum errors and their directions can be obtained from polar plots of the 
complex errors if desired. 


Characteristic Disturbances for Nonsymmetric Spacecraft 

Most disturbances acting on rotating spacecraft may be approximated by impulsive 
torques, step torques, step products of inertia, or variable products of inertia. For 
example, docking impacts and attitude control moments can be represented by impulsive 
or step torques, whereas crew or cargo motions result in either step or variable products 
of inertia. Other externally applied torques (such as the sinusoidal gravity gradient 
moments) are dependent on the particular spacecraft and orbital characteristics and can- 
not be defined without selecting a specific vehicle and orientation. 
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The effects of characteristic disturbing functions on the spacecraft motion are pre- 
sented in this section. Time solutions for the Euler angle and body rate errors are 
developed from the solution functions of tables I to V for arbitrary constant moments of 
inertia, and upper bounds for these variables are given. 

Impulsive torques .- Docking impulses caused by resupply and rendezvous vehicles 
or micrometeorite hits may result in impulsive torques acting on the spacecraft. 

Time histories: The corresponding time-dependent torques can be written as 

M = M x + iM y = (T x + iT y j 6(t) (47) 

and the corresponding forcing function is 


F= f 

■be 


T x 6(t) - XyTy6(t) 


+ Iy 


T y 6(t) + X x T x 6(t) 


(48) 


For arbitrary initial conditions, the total rate error may be found from equations (27), 
(42), and (48) by using the solution functions given in table I. The resulting expression is 


£2xy - &xy + 


'Vv 

T x cos Xt ^ sin Xt 

- X 


IvL 


T 

Ty cos Xt + sin Xt 


(49) 


where fijjy is given by equations (39) and (40). 

The attitude error can be similarly determined as 


a = a + 


- i 


Ty(cos Xt - cos at) + T x (--^ sin ^ + sin 
T x (cos Xt - cos at) - T y sin Xt + sin atj 


with 2 determined from equations (39) and (41). 


(50) 


Total errors: Conversion of the total angular and rate errors to exponential form 
leads to the complex vector representation 


fixy = Ojj + if^y = f^xy + B^e 


(51) 
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where 


^xy — + ^2© 


(52) 


and the complex constants B-j are 




“i- 1 


1 ; + 1 


lx K ly ^y 


B 2= f 


lx ^ ly \Z^y 


[(Wiofi) + ^y*fyo|/*y) >e 


Ix^xo^) + (*ytyofiyJ 


itan" 


VVo jfiy 


-it an 


Ix<W^ 


V^yo l/V 


■'■x^xol/^x 


> (53) 


r 


- : itan ^ 

T yl/*y 

B 3=<; 

V- 

if 1 1 

! |lx\pS ly^dl 

1/ ,— \ 2 t \ 2 1 

|( T x^j + ( T y/v) ) e 

Jxl/V 


b 4 = 


tyfiyl 


-itan 


-II 


T x^x) + ( T y^) 


T y^y 


[jx^Sc. 


The angular error oj can be developed directly by substitution of equations (52) 
and (53) into equation (18). The result is 


Q5[ = (<p + i0)e lcr * = Sj + i 


C 4 + Cge 1 ^* + C-e 1 ^-^ 


(54) 


where 


“i 


= 5e* a * = Oq + i 


Cl + Cae 1 ^^ 4 + C 3 e 1 ^ cr- ^ t 


(55) 
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and the constants are 



Initial error contribution: The total errors fix y an< ^ which correspond to 

the initial conditions fix 0 , f2y 0 , a, cp 0 , and 0 O are considered first. A simple geo- 
metrical interpretation of these error traces is possible. For fixy this interpretation 
follows from the trace of the velocity error in the XY body axis plane, as shown in 
figure 5. The path described by the tip of the Oxy vector is an ellipse in the body- 
fixed plane. 
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The characteristics of this ellipse 
are derived from an examination of equa- 
tions (52) and (53). The semiaxes a 
and b are determined as 


a = ^=^Ix°xo^x) + ( I y^V°'/v) 


) (57) 


I 2 2 

b = i + ( I y f2 y°'^) 


Y 



and the angular position of the rate error 
vector is given by 


Figure 5.- Rate error trace for initial conditions. 


$ a = tan V JL1 l 2E tan 


Uyfiy 


\lx^xol/^xy 


(58) 


The quadrant for the angles corresponding to the inverse trigonometric functions 
tan - *^^ in equation (58) and in all subsequent equations is determined by the sign of 
the numerator and denominator of the term in the braces. When both numerator and 
denominator are positive, the angle is in the first quadrant; when the numerator is posi- 
tive and the denominator is negative, the angle falls in the second quadrant; when both 
numerator and denominator are negative, the angle is in the third quadrant; and when the 
numerator is negative and the denominator is positive, the angle falls in the fourth 
quadrant. 


The position of the major axis of the ellipse is determined by the relative magni- 
tudes of I x , and I z . If I z is a maximum inertia, then the maximum angular rate 
will occur about the axis corresponding to the larger of the inertias I x and Iy ; if I z 
is a minimum inertia, then the maximum angular rate will occur about the axis corre- 
sponding to the smaller of the inertias I x and Iy. The period of revolution is 
and the S2xy vector rotates in the direction of the precession rate X. When I z is a 
maximum inertia, this rotation is in the direction of spin; when I z is a minimum 
inertia, the rotation is against the direction of spin. 


2n 

X 
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The trace of the rate error vector can be directly compared with the results of 
Poinsot's geometric construction (ref. 12), in which the path of the instantaneous rate 
vector on the ellipsoid of inertia is called the polhode. For the present solution the rate 
vector is restrained to move in a plane normal to the Z-axis, which is a principal axis 
of the inertia ellipsoid. The polhode projection onto this plane has been developed by 
Thomson (ref. 6, p. 124) and yields a curve whose shape is defined by the relation 

^x(lx^x) 2 + ^-y^y^y) 2 = ^x(lx^xo) + ^y(ly^yo) 2 (59) 

This relation describes an ellipse, with semiaxes given by equations (57). Since the 
polhode projection is proportional to the rate vector trace derived in this analysis, the 
approximate solution exactly represents the spacecraft rates when the variation in the 
spin rate is negligible. 

The angular trace with respect to the Xj- and Yj-axes is illustrated in figure 6. 

The path described by the tip of the S[ vector is generated by a point moving on a 
displaced ellipse which, in turn, is rotating at the spin rate. 




Figure 6.- Attitude error trace for initial conditions. 



From equations (55) and (56), the center of the moving ellipse is located by the 
vector sum of the initial attitude error oiq and the initial angular momentum ratio term 

i x ^ xo + . The radius a shown in figure 6 is thus 

CTl z J 


a 




(60) 


while the semiaxes of the rotating ellipse become 


b - ““in ^I x fi xo j + (VVojpy) 

C = ai 1 + ( I y^yo\fiy'j 

y 


(61) 


and the precession of the attitude error vector within the ellipse is specified by the angle 



Xt + tan 


-lfh^yofty} 


(62) 


When X is rational, the path of the attitude error trace is closed and has a period of 
2k7r, where k is the least common denominator of a and X. 

The trace of the attitude error vector is in agreement with the general properties 
predicted by MacMillan (ref. 12) for the torque-free motion of a rigid body with respect 
to a unit reference sphere. This sphere was drawn about the fixed point of the spinning 
body as a center, and the motion of the body Z-axis about the fixed- momentum axis was 
then described by the trace of the Z-axis on the unit sphere. The vector trace, intro- 
duced in the present analysis, can be considered to be the projection of this Z-axis trace 
onto a plane perpendicular to the Zj-axis. 

It should be apparent that the Zj-axis, which is arbitrarily defined as the fixed 
space axis corresponding to the initial position of the Z-axis, does not generally coincide 
with the fixed- momentum axis. By assumption, however, the angle between these two 
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axes is small. Hence, the shape of the traces about the fixed-momentum axis should be 
approximately retained in the plane normal to the Zj-axis. The fixed- momentum axis 
appears as a displaced point on this attitude error plane. 

In figure 6, the fixed- momentum axis projects as the center, of the rotating ellipse. 
The attitude error oscillates between two concentric circles drawn about the ellipse 
center. The radii of these circles are given by the minor and major semiaxes of the 
ellipse. The similarity of this motion to that depicted in figure 61 of MacMillan's text 
is obvious. 


Upper bounds of the values for the rate and attitude error magnitudes, as developed 
from equations (45) and (46), are 


|°xy| 


lim 


T 

Ix Iy fiy 


Ix^ ^yfiy 


^ (lAofc) + (^y^yofiy) 


(63) 


and 


Hlim " ^o 2 + 0 o 2 + ^\Jfh^xof + (Iy%o) 2 + ^ 


4 —i 

V x y v x x 


_1 1 _ 

fiy 


&XO fix') + (ly^yofiy) 


(64) 


The rate limit (eq. (63)) gives the major semiaxis of the rate error trace and is equal to 
the maximum rate error. The attitude limit (eq. (64)) corresponds to the sum of the 
center radii and the major semiaxis of the attitude ellipse and is greater than or equal 
to the maximum attitude error. 


Several interesting trends may be observed from the geometrical development and 
the relations for the upper limits of the errors. When the spacecraft inertia I x (or Iy) 
approaches I z while the second inertia Iy (or I x ) remains different from I z , then 
the rate and attitude ellipses become very elongated. Small rate errors induced about 
the second inertia axis by impulsive torques or other disturbances can thus lead to large 
total attitude and rate errors. An example is a cylindrical configuration spinning about 
an axis normal to the axis of symmetry. 
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When the spacecraft inertia I x (or Iy) is very much larger than I z , excessive 
attitude errors are produced by small body rates and tumbling may occur. This result, 
however, is not surprising since the inplane angular momentum is now much larger than 
the spin momentum. Examples here are slender cylindrical satellites and missiles 
spinning about a minimum axis of inertia. 

It may be noted that the smallest errors are produced when both I x and Iy are 
much smaller than I z and the spacecraft configuration approaches that of a disk. 

The contributions of the errors Oxy and 5 to the limiting errors | ^xy|jj m 
and |a|n m for a given disturbance are omitted in the remainder of the analysis to 
avoid undue complications of the limiting error relations. These error terms could, 
however, be readily included if this is desirable for a particular disturbance. 

Impulsive torque contribution: If the initial error terms are taken to be zero, the 
total errors for the impulsive torques are equivalent to those for the initial rate error 
terms. The geometrical representation and the maximum error values for the vectors 
corresponding to the initial errors thus hold for the impulsive torques if f2x 0 is replaced 
b y T x/ J x> Oyo is replaced by Ty/ly, and <p Q and 6 Q are taken to be zero in equa- 
tions (57) to (64) and in figures 5 and 6. 


Step torques .- The spacecraft attitude control system and external sources, such 
as gravity gradients, may also exert torques about the body axes. 

Time histories: For this example, consideration is given to constant step torques 

of the form 

M = M x + iM y = (T x + iTy)u(t) (65) 

and the associated forcing function 


F = f 


T x 5(t) - AyTyUOoJ + j^jTy6(t) + X x T x U(t)j 


( 66 ) 


The body rate error is found from equations (21), (42), and (66) by substitution of 
the solution functions of table I, with the result 
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°xy = ^xy + ^\r- 


T X 

T x sin Xt X^"(l " cos **) 


+ iy 


V* 


T y sin Xt + — 


(1 - cos Xt) 


(67) 


The attitude error a determined in a like manner is 


a = a + 


ctIz 


—\^r- - cos Xt| - icos at 

Xyliy j <y 


+ T, 


— sin Xt - i sin at 


+ ION 


l4-(-^ - cos Xt] - 4 cos atl 


Mix 


+ T, 


- sin Xt - 4 sin at 
X o 


( 68 ) 


The initial error contributions f2xy and 2 are defined as before. 

Total errors: A transformation of the errors to polar form yields 


Oxy = Oxy + B 3 + B 4 e lAt + B 5 e“ lXt 


(69) 


with 




-itan 


-1 


Tx X X*X 




B 4= f 


B 5= f 


1 1 1 

,kv 4 

_Ix^ lyfiy_ 

1 Xy X x 

J 


-itan" 


T x^" 


Ty^j 


1 1 l 

MVl 

r 

*xv 

Jy/ 

h&k bfi y 

\j Xy a. x 

r 



J 


(70) 


and, similarly, 


~ ~ n „i(CH-X)t n _i(cr-X)t . n 

Qjj. — ofj- + Cjj + Cge + C^e + CgC 


(71) 
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with 


C5= {^ 


itan 


-1 


T x 2 + Ty 2 W 


X 


c 6 = \- 


c 7 = <- 


C8 = 


1 

1 

2I z ct 

fy 


/t 2 t ^ I 

P- + -X 

I Xy X x 


itan 


-1 


Zzfi 


TxV^i 




ifq n ^ 

Ty^ 

l 

r 1 ii 



2I z a 

& \I X X 

j Xy X X J 



vv* 



itan 


-1 


T y x y I y 


^x x x^x 


(72) 


The trace of the total rate vector due only to the step torques is shown in figure 7. This 
trace is now an offset ellipse which intersects the origin at time zero. 


The elements of the rate ellipse can 
be determined from equations (69) and (70). 
The center of the ellipse is located by the 
radius 


a = 


vVy i 


\ ^xlx j 


(73) 


and the angular coordinate 


<& a = tan' 


■ i( T y¥y 

bx X X^X / 


(74) 


Y 



Figure 7.- Rate error trace for step torques. 
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The semiaxes of the ellipse are found from 


b = 


c = 


Ix^ 


|T X * T y" 


1 1 

|rp 2 rp 2 

1 X i y 

IyV 


| Xy X x 

3 J 


( 75 ) 


and the angular location of the rate error vector is 


1 

tan V tan 

it - 

r 

\ T y\^y/jJ 


(76) 


The position of the major semiaxis of the ellipse is again dependent on the relative mag- 
nitudes of I x , Iy, and I z . The maximum semiaxis is parallel to the larger (if Iz is a 

maximum inertia) or to the smaller (if I z is a minimum inertia) of the inertias Ix and 
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Iy. The motion of the S2 X y vector is in the direction of X and has a period of 

An extension of Poinsot's development (ref. 12) to the motion of a rigid body under 
step torques appears possible. The fixed reference point, with respect to which the 
polhode projection is generated, lies along the maximum angular momentum vector pos- 
sible for the body. The angular accelerations vanish for steady spin about this axis of 
maximum angular momentum. By referring to equations (2) and (3), it can be seen that 
the associated coordinates for the fixed point are proportional to 


^xo — 


*x*x 


(77) 


fyo “ \ T 

y Xyly 


for the constant step torques. 

The shape of the polhode projection corresponding to this fixed point is defined by 
equation (59). Substitution of equations (77) into equation (59) yields the curve defined by 


^xfjx^x) 2 + Xy([yfiy) 2 - 


(78) 
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An inspection of figure 7 shows that the polhode projection is indeed represented by the 
rate vector trace. The fixed point coincides with the center of the trace ellipse and the 
equation of the ellipse becomes equation (78). 

The inertial attitude error corresponding to the step torques yields the trace shown 
in figure 8. This trace results from a point moving along an ellipse, which remains 
fixed with respect to the rotating radius of a stationary displaced circle. The motion 
begins at the origin at time zero. 



Figure 8.- Attitude error trace for step torques. 

The center of the stationary circle is located by the radial coordinate 


o*l V 


IT* ^ 4. T ^ 
A x + l y 


(79) 


corresponding to the ratio of the torque to the spin vis viva. The radius of the stationary 
circle is 

, 1/ T-_ \ 2 / T„ \ 2 

(80) 



and the center of the moving ellipse has the angular position 


<t>k = at + tan 


-l &ffi ). tan-lfe) 

V^x^x^x/ \Tx/ 


(81) 
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The semiaxes of the ellipse are computed from the parameters 


c = 


d = 


0l Z )JX y 


viz ife 

and the semiaxis |c| makes the angle 


T 1 2 T> 2 


rtf tf 

XxT Xv 


^ c =l-tan- 1 ( ryIyXy 
2 \TxIx^x/ 


( 82 ) 


(83) 


with the radius b of the stationary circle. The angular coordinate of the tip of the atti- 
tude error vector within the ellipse 



completes the development of a point on the error trace curve. 


(84) 


The motion is a closed curve when X and a are both rational. The corre- 
sponding period is given by 2kir, where k is the least common denominator of <x 
and X. 


The general properties of the motion can be readily interpreted from figure 8, if it 
is recalled that this figure represents the projection of the Z-axis trace onto a plane nor- 
mal to the Zj-axis. The fixed space axis corresponding to the axis of maximum angular 
momentum projects as the center of the stationary circle. The motion of the body is 
bounded by two circles, concentric with the stationary circle and tangent to the moving 
ellipse in the figure. The outer circle represents the maximum Z-axis excursion rela- 
tive to the fixed momentum vector and can only be approached from the inside. The 
inner circle represents the minimum Z-axis excursion and is approached from the out- 
side when the ellipse does not contain the fixed momentum reference point. When the 
ellipse contains the momentum reference point, then the inner circle is crossed by the 
Z-axis trace and is approached from the inside. 


Upper bounds of the error magnitudes are found from equations (45) and (46), with 
the result 
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These upper bounds consist of the sums of the radial vector magnitudes and the semi- 
major axis of the error ellipse and are obviously greater than or equal to the maximum 
error values. 

For cylindrical spacecraft spinning about an axis normal to the cylinder axis, large 
rate and attitude errors will be produced by torques applied about the cylinder axis. 
Conversely, torques applied about the normal inertia axis in the spin plane will have 
little effect on the spacecraft motion. 

Near-cylindrical spacecraft spinning about a minimum inertia axis, so that I z is 
much less than I x (or Iy), will now be stable if the applied torques do not approach 
the spin vis viva term I z a 2 and if no energy is dissipated. 

Step products of inertia .- Crew or cargo movements within the spacecraft may be 
represented by equivalent masses mj with variable position coordinates xj, yj, 
and zj. The movements of the equivalent masses fall into two categories, namely, non- 
periodic motions and periodic motions. The first category includes arbitrary non- 
periodic motions along linear paths to some final position. From previous results for 
symmetrical spacecraft, it appears that the largest rate and attitude errors for such a 
motion are less than or equal to those for instantaneous motion to the final position. The 
introduction of step products of inertia corresponding to the final position coordinates of 
the moving masses gives a limiting case for this type of motion. 

Time histories: The coordinates of the jth mass may thus be written as 


x j “ x jo 

y j ■ y j° > 
Zj = z JO U(t) 
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and the corresponding forcing function is 

F = - ilyzj^(t) + 0^6(t)J + Xyl xz j^(t) + °^U(t)J, 


iX z [e(t) + O^t)] - X x I yz [6(t) + o2u(tl 


( 88 ) 


where the products of inertia now take on the constant values 
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The solution for the rate error becomes 
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and the attitude relation yields 
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Total errors: The vectorial representation of the total errors reduces to 
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The trace of the rate error vector is illustrated in figure 9. Once more a dis- 
placed ellipse is obtained. However, the ellipse now does not intersect the origin. 



Y 



Figure 9.- Rate error trace for step inertia products. 
The radius to the center of this ellipse is 


(w)’ 

and the ellipse characteristics are found from 


-\ 



The major semiaxis of the ellipse is parallel to the larger, for I z > 1^, or to the smaller, 
for I z < Ik, of the inertias Ik- However, the body axis with the largest rate error 
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is determined primarily by the location of the center of the trace 
the counterclockwise motion is again the precession period 
motion is in the direction of the precession rate X. 


ellipse. The period of 
and the direction of 


To correlate figure 9 with Poinsot’s development, note that the coordinates of the 
fixed reference point for the polhode are proportional to 





(99) 


from equations (2) and (3). This shift of the reference point from the coordinate origin 
to the principal body axis indicates that the maximum possible angular momentum vector 
lies along the new extreme inertia axis. Steady spin, for which the polhode reduces to a 
point, is thus possible only about the new maximum or minimum principal axis of inertia. 
In figure 9, this principal axis passes through the center of the rate ellipse, as specified 
by equation (96). 


The interpretation of the polhode projection about the fixed reference point becomes 
somewhat more difficult. Two terms now contribute to the polhode, namely, the rotation 
of the extreme principal axis and the effective acceleration torque produced by introduc- 
tion of the step product of inertia. The rotation of the extreme principal axis yields the 
initial rates given in equations (99). The step introduction of the product of inertia 
yields the additional rate terms 
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The polhode projection (eq. (59)) corresponds to the rate vector trace in figure 9, if the 
reference body rates Qxq and Qy 0 are given by the sum of equations (99) and (100). 

The angular position trace for the step products of inertia is presented in fig- 
ure 10. The trace is produced by a point on an ellipse, which remains fixed with respect 
to the rotating radius of a stationary circle centered at the origin of the inertial axis 
system. 
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Figure 10.- Attitude error trace for step inertia products. 


The radius of the stationary circle is 




and the ellipse center is defined by the angle 
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The ellipse semiaxes are found from 
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and the angle between the stationary circle radius a and the ellipse semiaxis | b | 
becomes 
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The position angle of the attitude error vector tip is now 


tan 

Xt-tan-lf Ixz ^|| 


\Iyz\[\x/ j 


with respect to the ellipse semiaxis |b|. 


(105) 


The period of the precessional motion is 2ki t, with k taken to be the least com- 
mon denominator of a and X for rational a and X. 


The general properties of motion for the step inertia products are similar to those 
for the step torques. However, the fixed space axis corresponding to the maximum 
angular momentum vector projects as the origin of the XjYj- plane, and the direction of 
the total angular momentum vector for the spacecraft is not changed in inertial space 
during the step crew motion. The spacecraft motion is bounded by two circles, drawn 
with center at the origin and tangent to the moving ellipse. The nature of these bounda- 
ries has been discussed in the step torque analysis. 

Upper bounds of the errors are 


| °xy| lim 
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for the rate vector and 
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for the attitude error vector. 
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As before, cylindrical spacecraft spinning about a normal axis lead to large errors 
for small products of inertia in the plane corresponding to the two large inertias. Some 
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differences in the response for the step inertia products and that for step torques can, 
however, arise after the removal of the disturbance. For the step torques this removal 
can occur when the rate error vector passes through the origin of the body axis system, 
so that the only residual error is a constant attitude error corresponding to the attitude 
at the time of disturbance removal. This fact may be of use in the design of pure atti- 
tude control systems for spacecraft which use constant torque pulses to reorient the 
spacecraft. 

For step products of inertia, removal of the disturbance could also null the rate 
errors if done when the rate error vector passes through its initial position. In practice, 
the determination of this position does not appear feasible without very exact values for 
the spacecraft and disturbance characteristics. The elimination of the body rate errors 
by the timely removal of the product of inertia is a very complex task and, in general, 
leads to both residual rate and attitude errors. 

Variable p r odu c ts of inertia .- The second category of mass movements within the 
spacecraft involves periodic motions, such as mass transfer along a circumferential 
path. The uncontrolled motion of the spacecraft is now similar to that of a spring-mass 
system with a periodic forcing function. The amplitudes of the rate and attitude errors 
are correspondingly multiplied by a magnification factor and resonance may occur for 
particular mass transfer rates. The determination of the effects of periodic mass 
motions within the spacecraft is thus an essential prerequisite to the analysis and selec- 
tion of the spacecraft control system. 

There are, unfortunately, an infinite number of possible periodic crew motions. 

The best approach to a study of these motions may be the formulation of a general 
forcing function, which permits the development of stability criteria and rate and attitude 
errors for a number of representative periodic motions. Particular time histories can 
then be developed for special cases of the general forcing function. 

To arrive at such a general forcing function, assume that the periodic mass motions 
involve transfer of a single mass and take place in a spin plane perpendicular to the 
Z-axis or along a line parallel to the Z-axis. If a sinusoidal variation is adopted for the 
associated mass position coordinates, the complex forcing function takes the form 

F = Eq cos pt + sin pt + E 26 (t) + Eg 6 (t) + E 4 &(t) (108) 

The coefficients Ej are complex constants which must be determined for par- 
ticular mass motions, and p is the frequency of the motions. 

A number of characteristic motions included in equation (108) are listed in table IV 
along with the corresponding coefficients Ej. As seen from the table, the forcing func- 
tion equation (108) comprises linear periodic motions parallel to each spacecraft axis 
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in the spacecraft reference planes XZ and YZ and circumferential motion in a plane 
perpendicular to the spacecraft Z-axis. Other mass motions can be constructed by com- 
bining the forcing functions in the table and by adding the forcing function equation (88) 
for the static products of inertia with appropriate values of I xz and Iy Z to the results. 
Any linear oscillation in a spin plane perpendicular to the Z-axis or along a line parallel 
to the Z-axis can be developed by this method. 

Only the motions described in table IV are considered as examples for the present 
analysis; the results obtained for these motions can be readily applied to more complex 
motions obtained by linear combinations of the forcing functions given by the table. 

Stability trends: An assessment of the stability trends for the variable products 
of inertia can be obtained by examining the forcing function terms in equation (108) and 
the solution functions of tables I and III. The uncontrolled spacecraft motion is unstable 
when the roots of the governing equations (29) and (34) contain real, positive, nonzero 
terms or when the forcing functions produce resonance conditions. From table I, it 
is apparent that the solution functions for 5(t), 6(t), and '5(t) contain only constant 

and periodic terms. These terms thus cannot cause divergence of the rate and attitude 
errors. 


The solution functions for cos pt and sin pt can, however, contain divergent 
terms and may lead to continuously increasing errors for special frequencies of the 
periodic motions. These special frequencies are |p| = |a| and |p| = | A | . ' When 
| p | = | cr J, the rate errors remain bounded for all finite values of Eq and Ej, but 
precession of the spacecraft may result unless 


Eg + iEj = 0 
Eq - iEj = 0 


(for p = o)~l 

> (109) 

(for p = -o)j 


Since equations (109) hold for all the values of Eq and E^ in table IV, this precession 
does not occur for the examples considered here. 

When |p| = |x|, both the rate and attitude errors tend to diverge for nonzero values 
of Eq and Ej. Mass motions with this period thus exhibit definite trends toward 
instability and should be avoided. 

There is, of course, one other instability that may occur for the present solution. 
From the governing equation (29) for the body rates, it can be shown that the rate error 
diverges when X^ < q ? so that I x > I z > Iy or Iy > I z > I x . This condition results when 
the Z-axis is an intermediate axis of inertia and agrees with the well-known fact that the 
undamped spacecraft spin is stable only if the spin axis is an axis of maximum or mini- 
mum inertia. In terms of the moving mass parameters, 
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may be taken to be the alternative conditions for instability corresponding to the assump- 
tions of this analysis. 

The conditions in this section serve to indicate possible instabilities for the 
rotating spacecraft. Since the rate and attitude errors for these instabilities rapidly 
exceed the small angle and rate assumptions, time histories for these motions are not 
discussed. If the unstable motions do occur, exact computer solutions should be used to 
assess their effects. 

Time histories: From equation (108) and the solution functions of tables I and in, 

the rate error becomes 


^xy ~ ^xy + 



cos pt + 



sin pt + E 4 6(t) + E 5 cos At + Eg sin At 


where 


E 5 = e 3 



E 6 = i 


(e 2 - x 2 e 4 ) - 


"A 





( 112 ) 


( 113 ) 


46 



and the attitude error is 
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Total errors: The vectorial representation of these errors yields 
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where 
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b 7 

C 8 = -a^p 

The total error relations given by equations (112) to (119) may now be evaluated for 
particular motions by substitution of the corresponding coefficients Eq, Ej, E 2, E3, 

and E4 from table IV. Only the vectorial representations and the upper bounds of the 
errors will usually have to be considered in an assessment of the effects of mass motion 
on the spacecraft motion. 

To determine the vectorial traces, values for E4, B4, B5, Bg, and B7 are 

needed. These expressions have been developed and are shown in tables IV and V for 
the motions described in table IV. The polar plots of the rate errors now follow directly 
from equation (116) and are presented in 
figure 11 for a number of typical motions. 

These traces are generated by a 
point on an ellipse, whose center moves 
along a second ellipse centered at the 
origin of the body axis system. Figure 12 
shows this development of the rate vector 
trace. 

The direction of the ellipse semiaxes 
| a | and |c| is first determined by the 
relation 
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Rate error trace 
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Figure 12.- Rate error trace development for periodic 
inertia products. 
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The position angles for the fixed and moving ellipses then become 



(121) 


The fixed ellipse has semiaxes |a| and |b|, and the moving ellipse has semiaxes |c| 
and |d|. The constants a, b, c, and d as determined from equations (116) and (117) 
are 

a = Bg + B 7 I 


b = Bg - B 7 


c = B 4 + Bg 




( 122 ) 


d = B 4 - BgJ 


where the values of the Bj terms are taken from table V. 

The direction of motion is determined by the signs of p and X, and the rate vec- 
tor describes a closed curve when p and X are both rational. The period of motion 
is 2kn, where k is the least common denominator of p and X. 


The rate error trace results 
from the oscillation of the mass with 
respect to the geometric body axes 
(p-ellipse) added to the precession 
of the geometric body axes about the 
principal body axes (X- ellipse). The 
fixed reference point for the trace 
falls on the origin of the geometric 
body axis system since the space- 
craft momentum must be conserved. 

Characteristic traces of the 
attitude errors for the periodic 
motions are illustrated in figure 13. 
These traces are produced by two 
ellipses which rotate at the spin rate, 
as sketched in figure 14. 



Figure 14.- Attitude error trace development for periodic 
inertia products. 
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The center of rotation is defined by the coordinates 
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The direction of the ellipse semiaxes |b| and |d| is determined by the angle at, 

measured from the positive Yj-axis when <3> a = 0 or 77 and measured from the negative 

Xj-axis when <$ a = — or -s-. The angular coordinates within the ellipses now become 
1 2 ^ 


= tan 


-1 


^ tan pt 


■§ tan Xt 
d 




(125) 


<£<3 = tan“ 

and the ellipse semimajor axes |b|, |c|, |d|, and |e j are derived from 
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The motion is a closed curve when o, p, and X are all rational. The precession 
period is 2k7r, where k is the least common denominator of a, p, and X. 

The attitude trace for the spacecraft with periodic inertia products exhibits a 
change in the direction of the angular momentum vector in inertial space. This rota- 
tion of the momentum vector to the center of the p-ellipse results from the initial mass 
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acceleration terms, which exert a torque on the spacecraft. The motion, in body coor- 
dinates referred to the fixed momentum reference point, can be visualized as the sum of 
the mass oscillation with respect to the body axes (p-ellipse) and the precession of the 
body axes about the principal axes (X- ellipse). In the intermediate inertial coordinates, 
this motion is rotated through the angle crt. 

Values for the upper bounds of the rate and attitude errors may be calculated from 

|«xy|lim = | B 4| + |B 5 | + | B 6 1 + |B 7 | (127) 
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by using the coefficient equations for E4 and Bj from tables IV and V. 

From an examination of these coefficient equations, it can be seen that the maxi- 
mum error relations for small disturbance frequencies p yield the error limit terms 
corresponding to the introduction of step products of inertia. Both principal-axis- 
rotation and acceleration terms result for the circumferential motions; only principal- 
axis-rotation terms appear for the radial and vertical oscillations. A first estimate of 
the limiting errors for the periodic inertia products can accordingly be obtained from the 
appropriate step inertia product terms, when the disturbance frequency | p | is much 
less than the precession frequency |A|. 

As the disturbance frequency increases, the error limits also increase. As 
expected, divergence of the errors is predicted when |p| = |x|. For a further increase 
in the disturbance frequencies, the error limits continue to increase. 

When the disturbance frequency |p| is much greater than the spin rate o, the 
error relations for the periodic inertia products become directly proportional to |p|. 
Doubling of the disturbance frequency thus doubles the resultant error limits, and large 
errors can be introduced by small, rapidly oscillating masses which may occur in 
onboard motors, generators, or other equipment. 

Several additional trends are indicated by tables IV and V. For circumferential 
motion, maximum errors result when the sign of the angular velocity p of the motion 
coincides with the sign of the precession rate X. Motion at a negative spin rate 
(p = -0) eliminates all except the initial acceleration effects. Motion at a negative or 
positive spin rate (p = ±0) also nulls the errors caused by vertical mass oscillations, 
but does not significantly affect the errors caused by the radial mass oscillations. 
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Characteristic Disturbances for Symmetric Spacecraft 

A large number of spinning spacecraft will be symmetric about their spin axis, so 
that I = I x = Iy. These spacecraft include rockets and ballistic missiles which are spin 
stabilized to maintain their flight-path angle under initial body rates, unmanned satellites 
which are spin stabilized to maintain a fixed inertial position for communication and 
observation purposes, and large manned space stations which provide an artificial gravity 
field for their crew. A reduction of the general solutions developed in a previous section 
to the special case of symmetric spacecraft would accordingly have many applications. 

Most of the resultant solutions have been previously obtained by various approxi- 
mation and numerical integration techniques and are scattered through the literature 
(refs. 4 to 20). The results in this section thus make no claim to originality but do 
accomplish two important objectives. 

The first objective involves the determination of the form of the geometric error 
traces and of the maximum error limits for the various disturbances. These important 
properties of the motion have been only partially treated in the literature and tend to be 
obscured by the component form in which past solutions are primarily expressed. The 
simple trace geometry that results from the complex vector representation of the pres- 
ent analysis should be of considerable value. 

The second objective is the comparison of the approximate solution with applicable 
previous results. This comparison points out the principal differences between the pres- 
ent and past solutions and summarizes the trends of the motion with variations in the 
disturbance and spacecraft characteristics. 

The discussions of the polhode projections and of the motion representation by 
means of the unit sphere are not repeated here. If desired, these relations can be readily 
deduced from the analysis for the nonsymmetric case and the equations that follow. 

Impulsive torques .- For impulsive torque disturbances of the form 

M = M x + iM y = T6(t) = (T x + iT y )6(t) 
the forcing function (eq. (48)) becomes 


F = |[b(t) + iX6(t)] 

The total rate and attitude error equations (49) and (50) yield 


Oxy = £2xy + j{cos Xt + i sin Xt) 
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Initial error contributions in equations (130) and (131) are 

fixy = ^xyo( cos + i sin Xt) 


(132) 


and 


S = o^(cos at - i sin at) + i ^— ( cos _ 1 sin at) - (cos At + i sin At)j (133) 


from equations (40) and (41). 

The total error vectors for the impulsive torques reduce to 
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where the initial condition vectors are 
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The rate vector trace for sym- 
metrical spacecraft with initial rate 
and attitude errors is shown in 
figure 15. The curve traced out by 
the C&xy vector is a circle whose 
radius is the magnitude of the initial 
rate error vector £2xyo- The rate 
vector fixy rotates with the preces- 
sion rate X to generate the error 
envelope. 


The attitude error trace, shown 
in figure 16, is also a circle. The 
center of this circle is determined by 
Iflxyo 


the vector + i 


ctIz 


and the 


radius of the circle is the ratio of the 


Y 



Figure 15.- Rale error trace for initial conditions and 
symmetric spacecraft. 
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Figure 16.- Attitude error trace for initial conditions 
and symmetric spacecraft. 


2tt 

inplane and spin momentum. The period of motion is — , , . and the attitude error vec- 

+ A 

tor moves in the direction of spin. 

Maximum error values can be derived from figures 15 and 16. The result is 
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( 139 ) 


From equation (139), it can be seen that the spacecraft inertia ratio I/I z does not 
affect the maximum rate error. The inertia ratio does, however, enter into the relation 
for the maximum attitude error which decreases with a decrease in the inertia ratio. 


56 


Spacecraft whose inertia ratio approaches that of a flat disk thus yield the minim um 
total attitude error for a given set of initial error values. 

The error relations developed here have been partially described by Leon (ref. 5) 
for a spinning symmetric rocket, and the applicable present results agree with his con- 
clusions. His work does not, however, develop the detailed trace representation or the 
maximum error relations. Thomson (see p. 201 of ref. 6) includes the trace representa- 
tions for the initial errors, but his geometrical interpretation is incorrect. The initial 
attitude error c/q and the initial body rate fixy 0 are not generally orthogonal, as 
depicted in Thomson's work. 


The error traces and maximum error values for the impulsive torques can be 
T T 

again found by setting -p = ft x0 , -J- = Oy 0 , and <p Q = e 0 = 0 in equations (136) 
to (139) and figures 15 and 16. 

Step torques .- For constant step torques described by 

M = M x + iM y = TU(t) = (T x + iT y )u(t) 
the forcing function (eq. (66)) yields 



for the symmetrical spacecraft. 

From equations (67) and (68), the error time histories are 
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The rate error trace for step 
torques is given in figure 17. This trace 
is a displaced circle whose center is 
located by the vector The radius of 

the circle is — , and the period of the 
Al . ^ 

precessional motion is — . 

I A 

Curves for the attitude error trace 
are depicted in figure 18. This trace is 
generated by a point on a circle of radius 


whose center moves along 


the circumference of a displaced circle 

J t 2 fp 2 

with radius " ■ x y The center of 


this displaced circle is located by the 
T 

vector and the period of motion 

a 2l z 

is 2 k 7 t, where k is the least common 
denominator of a and X for rational 
cr and X. 


Y 



Figure 17.- Rate error trace for step torques and 
symmetric spacecraft. 
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The maximum values of the errors become 


|^xy| 


_ 2 V 


+ Tv 


max 


l|cr 


(145) 


and 


a 


lim 


I 7 + 1 


Ir. ~ I 


+ 1| 


frJ + Tv 


- 2 Iz 


(146) 


The error vector relations for the step torques and symmetric spacecraft agree 
with Thomson's results (ref. 6, pp. 198 and 207). Maximum attitude error limits for 
this disturbance have also been developed by Suddath and include the residual errors 
after removal of the disturbance (ref. 14, p. 8). His limits, which are smaller than or 
equal to the limits obtainable during the torque application, do not represent the worst 
case and are thus somewhat misleading. 

Note that the error limits are inversely proportional to the inertia difference 
|l z - 1 1 and predict minimum errors for a disklike configuration. When applied to a 
sphere the limits predict divergence, confirming the statement that a sphere cannot be 
spin stabilized. 


Step products of inertia.- When step products of inertia of the form 


A xz 


m j x jo z jo " m s x so z so 


11 

X 

j=l 

n 

lyz = ^ m jyjo z jo - m sy S o z so 

j=l 


(147) 


occur in symmetrical spacecraft, 


I rz - *xz + ilyz 


(148) 


may be introduced to get the forcing function 


from equation (88). 



(149) 
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The error solutions (eqs. (90) and (91)) give 

Ic2 - X 2 


ftxy - ^Xy + 


(“1^) ( g2 x - ) ( cos + i sin xt) - ^ + i6(t) 


and 


a = a + a - (a - A)(cos At + i sin At)j 


for the symmetric spacecraft. 
In vector form, 


(I rz 

&xy ~ ^xy + 1 ~ 


c2_-J^ e lXt . + 15(t) 


and 


°i. 


= 3[ + i^^e i<;t -((;-A)e i < a+X ) t ] 


(150) 


(151) 


(152) 


(153) 


The polar graph of the rate error is illustrated in figure 19. The error trace is a 

a 2 I r7 

The radius of the circle 


displaced circle with center determined by the vector - 


> T 


is 



and the precession period of the motion is 


277 

A 


Y 



Figure 19.- Rate error trace for step inertia products and symmetric spacecraft. 
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The attitude error for the step product of inertia is sketched in figure 20. The 
error trace now is produced by a point on a rotating circle with radius 
2 2 

xz + lyz which moves around the circumference of a circle with center at 


Mb 


the origin and radius QjjyI X z 2 + lyz 2 * The Period of motion is 2kn, where k is the 

least common denominator of a and X, and the trace curve does not close unless a 
and X are rational. 


Maximum error values become 

!^ x y|max 

for the rate vector and 


2 cr 2 - x 2 

XI 




for the attitude vector. 


. , 

2ct - A 

Hlim = 

XI 




(154) 


(155) 


The errors predicted by equations (152) to (155) do not agree with previous results 
(refs. 4, 5, and 6) for product of inertia disturbances. This might be expected, however, 
since the previous analyses have neglected the energy associated with the introduction of 
the product of inertia. 

The limiting errors (eqs. (154) and (155)) diverge for a spherical configuration. 
Attitude error bounds for a long slender cylinder (X — -a) are three times as large as the 
attitude error bounds for a flat disk (X — a). Rate error bounds for these limiting con- 
figurations are, however, equal. 


Y 


I 


Y 


I 



Figure 20.- Attitude error trace for step inertia products and symmetric spacecraft. 
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Variable products o f inertia .- Forcing functions, stability criteria, and solutions 
for variable periodic products of inertia and symmetrical spacecraft retain the form of 
equations (108) to (119). The symmetric spacecraft does, however, permit a simpler 
combination of the inertia products and yields a better understanding of the effects of the 
spacecraft parameters on the error traces. 

Consider a product of inertia of the form 

I rz = l xz + ilyz = Q z o( x o + Wo) ( 156 ) 

as produced by the motion of a single mass. Circumferential, radial, and vertical oscil- 
lations which yield this inertia product can now be examined. 


Circumferential mass motion: The circumferential mass motion begins at the 
point x Q = r 0 cos pt 0 , y Q = r 0 sin pt 0 , z Q and continues at a constant angular rate p 
around the perimeter of the spacecraft. The associated position coordinates may be 
written as 

x = r 0 cos p(t + t 0 )l 


y = r 0 sin p(t + t 0 ) 
z = z 0 U(t 0 ) 

r o = \/ x o 2 + yo 2 


(157) 


The solution coefficients are taken from table V, with the result 


B 4 


Irz (ct + p) 2 - (a 2 - p 2 ) P 

I L * - P + 2_j 


and 


b 5 


Be 


Irzp 

nr 

Irz( CT + P)^ 
(A - P)I 


B 7 = 0 




J 


P!rz 

~2 g - A 

A 

I 

[g 2 - A^J 



c- = - ^zRg + p) 2 - (a 2 - p2) P ~ 1 

5 I L (g + A)(A - p) 2(g+A)Jl 

- Plrz [ 

6 2(g - A) I 


= (g + P)I r z 
(A - p)I 


(158) 


(159) 
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The rate and attitude error vector equations (116) and (118) reduce to 


^xy = ^xy + Is |^P> 2 x: (x 2 -P 2 ) + |] e ‘« - E ,-U* . fjLL f e ‘P* + la(,| ( 160 ) 


and 


0i l = a l- — 


**rz J 

pp(2a - X)" 1 

_L 

(cr+ p) 2 - (x2 - p2) p 

i(a+X)t 

1 1 

1 

% 

1 

CO 
1 

T 

(cr+X)(X-p) 2(a+X) 

e 


P 

i(a-X)t _ 

It + p 

i(a+p)t\ 

J(a - X)j 

t; 

L x - p 

i 


(161) 


The rate vector trace, illustrated in figure 21, is produced by an ellipse, whose 
center travels along the circumference of a circle. The radius of this circle is 


(a+ p) 2 

^xz 2 + tyz 2 

__ X - p ^ 

I 


(162) 


and the ellipse center is determined by the angle pt, measured from the axis corre- 
sponding to the radial inertia product. The ellipse semiaxis |b| also lies along the 
radial inertia product axis, and the precession angle of the rate vector tip within the 
ellipse is 


where 


b 


c 

define the ellipse semiaxes. 


= tan" ^ tan" ^Xtj 


(g+ P) 2 
X - p 

(g + P) 2 
X - p 





(163) 


(164) 


The attitude error trace in inertial space, shown in figure 22, is derived by a 
rotating ellipse whose center translates along the circumference of a displaced circle. 


The center of the displaced circle is located by 

/p(2a - X)' 


a = 


*xz + lyz 


(165) 
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and the circle radius becomes 


b = 



(166) 


The center of the moving ellipse is specified by the angle (a + p)t and the semiaxis |c| 
rotates at the spin rate. 


and 


The angle to the attitude vector tip and the semiaxes of the ellipse are found from 


*c 



(167) 


■\ 


(q + p) 2 - (x 2 - p 2 ) pX 

,T 2 T 2 

yixz + -tyz 

X - p q - X 

(q + X)I 

(q + p) 2 - (x 2 - p2) pa 

2 t 2 

ylxz + J yz 

L X - p + 

(q + X)I 


(168) 


respectively. 

Limiting error values can be computed from the upper bound relations 

\/lxz 2 + Iyzjj (169) 

and 


|°xy| 


lim 


(q + p ) 2 - (x 2 - p 2 ) 

X - p 


+ 

P 

+ 

(ff+p) 2 


2 


X - p 



p(2q - X) 


(q + p) 2 - (* 2 - p 2 ) , p 

1 

p 

Mlim “ 

q 2 - X 2 

+ 

(q + X)(X - p) 2(q + X) 

T 

2(q - X) 


q + p 
x - p 




(170) 


These results indicate that the rate and attitude errors for the circumferential motion 
are larger than those for the static products of inertia, for which p vanishes. The 
limiting errors increase with an increase in p and tend to diverge as p approaches 
X. For values of p greater than X, the error continues to increase with an increase 
in p. 

For very small disturbance frequencies (|p| « q), the upper error bounds are 
given by 



2 o 2 - X 2 

\|lxz^ + V z 

X 

I 


(171) 
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and 


a 


lim 


2ct - A 

)jlxz + 

X 

i 


(172) 


As might be expected, the maximum errors now resemble the error bounds for the 
step products of inertia. The resultant bounds are essentially independent of the dis- 
turbance frequency. 

For very large disturbance frequencies (|p| » o), 


| 4 |p|^xz 2 + lyz 2 

l iix y|lim * I 

and 


(173) 



p(5a - 3A) 

/ixz + lyz 

o 2 - A 2 

I 


(174) 


so that both error bounds increase linearly with jpj. 

Relations (171) to (174) can be used to determine error limits and their variation 
with the disturbance characteristics, when the absolute disturbance frequency is much 
less than the absolute precession rate or when the absolute disturbance frequency is 
much greater than the spin rate. 


Radial mass oscillation: The radial mass motion begins at the point 0,0, z Q and is 

defined by the coordinates 

x = x Q sin pt 

y = y 0 sin pt > (175) 

z = z Q U(t) 


The solution coefficients, taken from table V, now become 

r- _-\ 


b 4 

b 5 

Be 

b 7 


ipl rz 

1 o 2 + 2aA + p2 

I 

5 X2 - p2 


ipfrz 

21 

i(a + p) 2 I rz 
2 (A - p)I 

-i(q - p)^I rz 

2 (A + p)I J 


(176) 
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and 


'iPlrz' 

1 q2 + 2qA + p2 

(q + A)l 

2 A 2 - p2 


„ _ip(2a-A)I rz 
4 ~ (a2-x2)l 


C 5 

C 6 = - 
c 7 = - 
c 8 = 


ipl 


rz 


2 (a - A)I 

i(ff + p)I rz 
2 (A - p)I 

i(g ~ P)I rz 
2 (A + p)I 


J 


Rate and attitude vectors are 




j _ 2 ( 0 ^ + 2aA + p 2 ) 
a 2 - p2 


e iAt + e -iAt 


(177) 


and 


(q + p ) 2 

P(A - p) 


ipt 


(q - p) 2 


p(A + p) 


-ipt\ 


(178) 


~ plrz 


f 2(2q - A)1 / 1 \ 

1 2(q2 + 2qA + p 2 ) 

l[_q 2 - A 2 J 

xO" j 


s i(q+A)t 


-(^ i<a - x)t -[p^^ +p,l + [i 


o- P 
P(A + P) 


a i(q~p)t\ 


(179) 


as determined from substitution of equations (176) and (177) into equations (116) and (118). 

The rate error trace, shown in figure 23, is generated by two ellipse envelopes, 
whose semiaxes | a | and | c | lie along a line normal to the radial product of inertia 
axis. 
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Figure 23.- Rate error trace for radial mass oscillation and symmetric spacecraft. 
The characteristics of these ellipse envelopes are developed from 



and the precession angles within the ellipses become 



(180) 


(181) 


The attitude error trace is also derived from two ellipses and takes the form 
illustrated in figure 24. 
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Figure 24.- Attitude error trace for radial mass oscillation and symmetric spacecraft. 


The ellipse semiaxes | b | and | d | initially coincide with the radial-inertia- 
product axis and precess at the spin rate. The center of the inner ellipse is given by 



and the semiaxes for the two rotating ellipses become 



from equation (179). The precession angles within the ellipses are 



(183) 


(184) 
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The upper bounds of the error magnitudes become 


Mlim " 


1 + 


j _ 2 (o 2 + 2oX + p 2 ) 


X 2 - p 2 



(a + p) 2 

i 

(ct - p) 2 

l?l/ 

+ 

p(X - p) 

T 

p(X + p) 



I 2 + I 2 
X xz + x yz 


21 


and 


In/I - 

2(2ct - X) 

1 

i 

1 2(o 2 + 2aX + p 2 ) 


1 

Mum “ 

a 2 - X 2 

+ a+ X 

X 2 - p 2 


a - X 


(185) 


a + p 

i 

a - p 

p il 

p(X - p) 

T 

p(X + p) 



T 2 x 2 

l xz + L yz 


21 


(186) 


These error limits again increase as |p| increases and diverge as |p| approaches |X|. 
When the disturbance frequency is very small (|p| « oj, the approximate upper 


bounds for the rate and attitude errors are given by 


n. 


xyllim 


o£ 

X 


II 2 + i 2 

l A xz + 1 yz 


(187) 


and 


Ml 


im 


n 2 + i 2 

i xz + yz 


(188) 


These upper bounds correspond to the maximum principal-axis rotation for the inertia 
product and do not contain the disturbance frequency. 

When the disturbance frequency is very large p| » aj, 


I SI 


2 |p|v I xz + l yz 


xyilim 


(189) 


and 


Mi 


im 


2p(2a - X) 

\jlxz + lyz 

a 2 - X 2 

I 


(190) 


The upper bounds defined by equations (189) and (190) increase linearly with p 
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A comparison of equations (171) to (174) with equations (187) to (190) shows that 
the errors for the radial mass oscillation are smaller than those predicted for the cir- 
cumferential mass motion, when the disturbance (p, I rz ) and spacecraft ( o , A, I) char- 
acteristics are equal. 


Vertical mass oscillation: The vertical mass oscillation comprises a periodic 
mass motion, which starts at the point x Q ,y o ,0 and is described by 


x = x Q 

y = y 0 > 

z = z Q sin pt. 


(191) 


The solution coefficients for this motion are 


b 4 = - 


ip(or2 _ p 2) Irz 

J? ~ P 2 )l 


and 


B 5 


= 0 


b 6 = 


b 7 = 


_ i(cr 2 - p 2 )l rz 
2 (A - p)I 
i(a 2 - p2)l r 

2(A + p)I 


^ _ -i P I rz 
° 4 (a + A)l 

_ iP^ - P 2 )l- 
5 (a + a) (a - p)i 


rz 


c 6 = 0 


c 7 = - 


i(tf - P)I rz 
2 (A - p)I 


c 8 = 


i(a + p)I. 
2(A + p)I 


rz 


J 


(192) 


(193) 
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The error vectors, as developed from equations (116) and (118), yield 


Oxy - Bay - - (X + p)e ipt + (X - p)e" lp ^j 




P_ 1 - ( *? ~ P 2 V (g-A)t] . l \ (v- P\ c i(o-+P)t _ ( o + £\ i(g-p)t 

+ x l\2 _ p 2/ e + 2 vx - py U + p r 


The rate error trace for the vertical periodic motions is given in figure 25. This 
trace is generated by a point on a circle, whose center traces out an ellipse in the 
XY-plane. At time zero, both the ellipse semiaxis | a | and the circle radius vector lie 
along an axis which is perpendicular to the mass position radius. 

The ellipse semiaxes are computed from 

-N 

x(cr 2 - p 2 ) V J xz + *yz 

A 2 - p2 I 


p(q 2 - p^^xz + *yz 


X 2 -p2 



Figure 25.- Rate error trace for vertical mass oscillation and symmetric spacecraft. 
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and the circle radius is 


c = 


pG* 2 - P 2 ) 

liP+W* 

X 2 - p 2 

HH 


The center of the circle along the ellipse is given by 


<h a = tan 


-li 


£ tan pt 


*] 


(197) 


(198) 


and the rate vector tip precesses at the rate X within the circle. 

The attitude error trace, sketched in figure 26, is produced by a circle which 
travels along the p- ellipse. Both circle and ellipse rotate at the spin rate a and are 
initially alined with the radial product of inertia axes. 

The center of the ellipse is determined from 


and the ellipse semiaxes are given by 


v 2 + hz 2 


(199) 


b = 


c = 



lx 2 r 2 

1 

b 

i'S 


A xz + *yz 

[x 2 - p 2 _ 

I 

ax - p 2 


I ^ + I ^ 
A xz + x yz 

X 2 - p 2 

I 


( 200 ) 



Figure 26.- Attitude error trace for vertical mass oscillation and symmetric spacecraft. 
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( 201 ) 


The center of the precessing circle is specified by 

<f>k = tan - 1 tan pt 
and the radius of the circle becomes 


' pf^-p 2 ) ' 


I ^ + I ^ 
A xz + A yz 

(q + X)(x 2 - p 2 ) 

I 


( 202 ) 


This radius moves with the precession rate X within the circular envelope. 

Values for the upper limits of the rate and attitude error magnitudes may be cal- 
culated by using 


and 


1 £2 



1 p 1 + i x - pi + lx + ph 

i 

| xy 

lim 

X 2 - p 2 

I I 2\l 1 1 1/ 



L xz 


+ I, 


•yz 


(203) 



a + p 
x + p 


+ 


q- p 


fxz 2 * hz 

X - p 

) 

I 


(204) 


The upper bounds of the errors increase as p increases, 
jpj = | X | is obvious from equations (203) and (204). 

If Jpj is small (|p| « aV these equations yield 


fix y lim 


q 2 ' 

L 2 T 2 

ybcz + A yz 

X 

I 


The instability trend at 


(205) 


and 



(206) 


so that both error bounds are approximately independent of | p | and contain only the 
terms associated with the principal- axis rotation. 

If | p | is large ^|p | » oj, the upper bounds reduce to 


£2 

*y lim 


2 P 




I 


(207) 
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and 


a 


[lim 



(208) 


as a first approximation. Both these error bounds are directly proportional to |p|. 

From equations (171) to (174), (187) to (190), and (205) to (208), it can be seen that 
the error bounds predicted for the vertical mass oscillation closely resemble those 
developed for the radial mass oscillation. The attitude errors for the periodic vertical 
motions are, however, somewhat lower than the attitude errors for the periodic radial 
motions when the disturbance and spacecraft characteristics are equal. 

Error trends: In concluding this discussion of the variable products of inertia for 
symmetric spacecraft, it should be recalled that the variation of the limiting errors 
with inertia distribution is similar to that described for the step products of inertia. 

Flat disk configurations yield the minimum errors, and spherical configurations are 
unstable. However, the magnitudes of the rate and attitude errors are now considerably 
larger than those of the errors produced by the step products of inertia. 


Controlled Spacecraft Characteristics 

Governing eq uations .- The motion of the controlled spacecraft can be defined by a 
method similar to that for the uncontrolled vehicle. The torques produced by the control 
system are now particular functions of the measured vehicle angular position and rate. 
These torques can thus be considered to be forcing functions applied to the uncontrolled 
vehicle equations. The solutions of the resultant differential equations yield the space- 
crafts angular position and rate errors, as before. 

The analysis begins with the selection of a control torque command and the develop- 
ment of the corresponding equations of motion. As an example, a linear control torque 
g is introduced as 

K lx fi x + K 2 x J °y dt + K 3 x 0 + • • • 


g “ g x J x + lg yV ~ 


+ i 


KlyS + K 2y 



dt + Kgy<P + . . . 


(209) 


where the are the physical control gains that must be provided by the stabilization 

system. The error signals are amplified by these gains, and a control moment, whose 
value is equal to the sum of the amplified error terms given by equation (209), is then 
applied to the spacecraft. 
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Particular nonlinear control torques, which lead to governing equations of the form 
discussed in the literature (refs. 21 to 24), could also be considered. The present 
application is, however, restricted to a discussion of the linear control functions in equa- 
tion (209). These linear control laws can be readily mechanized and allow a simple 
interpretation of the mechanics of motion for various types of sensor inputs (such as 
those derived directly from rate gyros, stable platforms, or Euler angle computers). 


In an analysis of the spacecraft’s stability, it is easier to deal with nondimensional 
control gains which may be defined by 

A 


u K lx 

lx XI x 


k iy = xf 


k„ 

2x AXyly 


_ K 2y 
xlx 


k 2y ~ xx j- 


w -%c 

3x axiv 


k_ _f?y 

3y aXI v 




= X 


so that g x and g y become 

k lx n x + k 2x(^) J °y dt + CTk 3x e + 

k lyf2y + k 2y(¥^y °x dt + ak 3y ^ + 


g y = X 


The equations of motion can now be written in the form 


hx + rf-^% = gx + f x 


% ■ - gy + f y 




( 210 ) 


( 211 ) 


(212) 


where 


n 

V 1 

^ m j( z j x j • x j z j) ' m s( z s x s - x s z s) 

a+ Z ■ y fi) - m s( z sy s - y s z s) 

i =1 

i-. j 


f x - — <M X + CT (ixz “ lyz°) + 


J-<|m x - ^ Z j (a 2 y ] - 2oxj + + z^ + m s z s (o 2 y s - 2ok s + y s ) + z s y s \ 

j= l J J 


( 213 ) 
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and 




n 

I m i( z 3^ 

j=i 



m s( z s y s - y s z s)j 


a+ Z m i( X 3*3 

3=1 


s j X j) - m s( x 




, jfj(° f2x 3 + 2o s , j- 5 i)- 



m s ^(o 2 x s + 2 °y s - xs) 



(214) 


with the spin rate a taken to be constant and positive. 
The Euler angle relations may be expressed as 



(215) 


and the simultaneous solution of equations (212) and (215) specifies the spacecraft's 
motion. 


Control requirements .- In practice, the spacecraft's rate and inertial position 
errors must be kept within specified dead bands which are determined by the spacecraft 
mission requirements. Control of the spacecraft to these accuracies may be provided by 
a reaction jet system and a momentum storage system, and the torques that must be pro- 
duced by the system actuators can be readily determined if the required control torques 
about the body axes are specified. With a reaction jet system, these torques are gener- 
ated by variable-mass-flow or pulse- modulated jets, and for a momentum storage system 
composed of control moment gyros or reaction wheels the corresponding torque com- 
ponents along the gyro gimbal or reaction wheel axis are computed and applied. Concern 
is given to the actual mechanization of such systems in a subsequent section; for the time 
being only the body-axis torques necessary to stabilize the spacecraft are developed. 

Since the magnitude of the angular error in inertial space is equal to the magnitude 
of the body-referred angular error, the damping of a and S2xy will assure the adequate 
stabilization of the spacecraft with respect to both body-fixed and inertial frames. The 
problem is thus reduced to the determination of control torques that will damp a and 
C X y to zero or to small steady state values. 

Control la w formulation .- To investigate particular control laws, the corresponding 
control torque functions g x and gy must first be specified. The associated governing 
equations follow from equations (212) and (215). The stability regions for the governing 
equations can then be defined by making use of the conditions developed by Routh (ref. 11) 
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and Hurwitz (ref. 25). If the selected control torques allow stable solutions for the 
Euler angles or body rates, time histories and complex error solutions can be found by 
the Laplace transform technique or by numerical integration of the linearized equations. 
If no stable solutions are possible, the control torques can be rejected immediately. 


To illustrate the applications of this technique, several example control laws are 
analyzed. 

Pure rate control law: For a pure rate control law, consider 




K lx 

g x - -j— - £2x ~ ^lx^x 

g y * ' ^ly^y 




so that the moment equations yield 


1 

^x " ^lx^x + = f X 

% “ ^ly^y _ f-^^X = % 


(216) 


(217) 


The Laplace formulation of the corresponding complex rate error is 


s 2 “ A ( k lx + k ly) s + a2 ( 1 + k lx k ly) fi xy( s ) = v ( s ) 


where the transform of the effective forcing function is given by 


(218) 


V(s) s V x (s) + iV y (s) = 


( S - Ak ly) fi X0 “ (^)%0 + (' 


+ 1 


(s - Ak lx )«y 0 + ^n X 0 + ( s - Ak ix) f y( s ) + 


- Ak lj^ f x (s) ” 




(219) 


The time solution of equation (218) consists of the sum of a general solution, for 
which f x (s) and fy(s) are set equal to zero in equation (219), and a particular solution 
of the complete equation (218), for which f x (s) and f y (s) are specified for the applied 
disturbances. The functions f x (t) and f y (t) are the explicit, continuous functions of 
time defined by equations (213) and (214) and do not contain the rate or attitude errors. 
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The particular solution of equation (218) is directly dependent on the disturbance 
under consideration, but characteristic trends for this solution can be indicated when the 
general solution is a damped vibration. Step functions in equations (213) and (214) lead 
to constant residual rate and attitude errors; impulse functions and their derivatives 
lead to damped transient rate and attitude errors which approach zero as time increases; 
and sinusoidal forcing functions produce residual sinusoidal rate and attitude errors. 

The amplitude of these residual errors is reduced with an increase in the damping ratio. 

The actual development of time histories for the various applied disturbances is 
not attempted here. The primary tasks of the control system are the minimization of 
errors during a disturbance and the elimination of residual errors after removal of the 
disturbance. Both these tasks can be accomplished by the selection of stable gains that 
yield large damping ratios consistent with realistic control systems and do not require 
the development of time solutions. 

Stable gains are gains for which all roots in the general solution have negative real 
parts. The characteristic equation for this general solution is 

s 2 - X(k lx + kjyjs + X 2 (l + kj x ki y j = 0 (220) 

and stability of the complex rate error requires that 

-N 

-Xk^ > Xk^ 
k lx k ly > "1 J 

It should be noted that equation (220) has the form 

p o 

s^ + 2rjja)^s + ajjq- = 0 

The damped natural frequency may be expressed as 

= + k lx k ly (223) 

and the damping ratio r-p and time constant t^ are, respectively, 


( 221 ) 


( 222 ) 


and 


D 


k lx + k ly 
2^1 + k lx k ly 


sgn X 


(224) 


t = 1 -2 

D r D w N X ( k lx + k ly) 


(225) 
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Selection of the control gains is generally based on a desired time constant and 
damping ratio. This selection makes use of the relations 


and 


where 



(226) 


(227) 


(228) 


for the specified real values of kj x and kjy The stability conditions of equa- 
tions (221) are automatically satisfied by equations (226) and (227) for positive real 
values of the damping ratio and time constant. Control gains for particular damping 
characteristics can thus be determined directly from equations (226) and (227). 

The resulting regions of stability for the control gain functions kj x sgn X and 
kjy sgn X are illustrated in figure 27(a). The gain functions must be in the area bounded 
above by the rectangular hyperbolas k^ x k^y = -1 and the straight line kj x = -k^y. Of 
particular interest is the fact that either of the two gains can be zero. This means that 
damping of the spacecraft's rates is possible with torques applied about a single space- 
craft axis and derived from a single rate gyro for that axis. Hence, a control system 
providing rate control torques about two axes has inherent redundancy in the event of 
failure of one of the system torquers. 

An indication of the damping characteristics obtained for the pure rate control is 
given in figure 27(b). This figure shows the overdamped, critically damped, and damped 
vibration regions corresponding to the stable control gains. Critical damping occurs 
when rjj = 1 and yields the conditions 

kfx = kj_y + 2 (229a) 

” k lx = k ly - 2 (229b) 

from equation (224). Stable gains that fall outside of the straight lines defined by these 
conditions will yield overdamped spacecraft motion; gains that fall between the straight 
lines will yield vibrational, damped spacecraft motion. 
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(d) Rate integral gains for small rate gains. 


(e) Attitude gains for small rate and attitude gains. 


Figure 27.- Stability characteristics for the controlled spacecraft 





From equations (224) and (225), it is apparent that the time constant is smallest 
when the two rate gains have the same sign. To optimize both the time constant and the 
damping ratio, one of the two gains can be selected to be zero. Single-axis rate control 
for the spacecraft should thus be quite efficient. 

The transformed complex position equation is 


(s + io) 


s 2 - \(k lx + k ly )s + X 2 (l + k lx k ly ) 


a(s) 


= V(s) + 


s - X^kj x + kj y ^ + A ^1 + kj x kj y j 


(230) 


Since this equation contains a purely imaginary root, the solutions for the complex posi- 
tion error are neutrally stable. The pure rate command is thus limited in its usefulness 
to those applications in which only rate damping is needed. An example of such an appli- 
cation would be a manned space station, where the functions of rate damping and attitude 
control are often provided by different subsystems and where the rate control law is 
used to command an onboard momentum storage system. 

From equation (218), the rate error may be written as 


S^xy : 


J(si+ s 2) 1 


X ( k lx ■ k ly)°xo - HH^y. 


sinhj 


( S 1 


•)t + 


^xo( s l “ s 2. 


coshfs^ - S2jt> + 


x(k ly - k lx )fiy 0 





( 231 ) 


where 


S l = 9 


( k lx + k ly) + )/( k lx “ k ly) " 




s 2 =| 


( k lx + k ly) - \|( k lx - k ly) 4 _ 


(232) 


The response of the spacecraft is analogous to that of a spring- mass system with forced 
vibration and damping. The free vibration term corresponding to the initial conditions 
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approaches zero as time increases, and the forced vibration term corresponding to the 
applied torque and mass motion effects is multiplied by a magnification factor or is 
damped to zero. The magnification factor is a function of the control gains and 
decreases as the terms kj x and kjy take on larger stable values. The attitude error 
for the pure rate control becomes 




+^1 

( s - k ly) f x( s ) - (^~) f y (s ) 


(s - Sl )(s 


+ i 


|( s - k lx) f y( s ) + 
2)( s + ix ) ~ 




f x( s )| 


( 233 ) 


where Sj and S 2 are given by equations (232). 

From inspection of equation (233), it is apparent that the initial errors (Xq and 
^xyo contribute both free and damped vibration terms to the spacecraft's motion. 

Since the damped terms vanish for large time, only the purely oscillatory terms need be 
considered. The amplitude of the motion due to the initial position error is not 
affected by the pure rate control law, as would be expected from the results of the sta- 
bility analysis. The amplitude of the motion resulting from initial rate errors can, how- 
ever, be reduced by proper choice of the control gains for the particular spacecraft 
under investigation. 

In summary, it may thus be concluded that the pure rate control law is adequate 
for the damping of the spacecraft's angular velocities. If no initial attitude errors or 
reorientation requirements exist for the spacecraft, it should also be possible to select 
control gains which will hold the spacecraft to small oscillations about its initial position 
in the presence of crew motions and other internal disturbances. This function is par- 
ticularly important in spacecraft with solar cell panels or similar equipment, which 
must be maintained approximately in a given inertial direction. The effectiveness of the 
attitude hold mode for the rate control law should, however, be checked by substitution 
of "worst case" forcing functions for the spacecraft into equation (233). 
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Rate plus rate integral control law: When the control torque is derived from both 
rate gyros and integrating rate gyros, 


g x = £( k 1x«x + k 2 x J % dt j = x k lx^x + k 2x(^j J % dt 

gy = + K2y ^ fijj d^ = X kjyfiy + 1^2y ^ Jy ^ ^ ^X 


(234) 


with the governing equations 


fix — ^^lX^X t j fly — XlC2x ^ % dt 

f2x + ^2y J"* ^x d ^ 


• / X-y'I'sA 

S2y - %y^y - K— ) 


= U 


= fv 


(235) 


The characteristic equation for the general solution for the rate error is given by 

s 4 “ X ( k lx + k ly) s3 + + k lx k ly) s2 + x3 ( k 2y " k 2x) s - X<lk 2x k 2y = 0 

and the stability conditions are 


(236) 


~\ 


-Xkjy > Xkj x 

k lx k ly > “* 

_k 2x k 2y > d 
Xk2y > 

( k lx + k ly) k 2x k 2y > ( k 2y “ k 2x)[( k 2y ~ k 2x) + ( k lx + k ly)(* + k lx k ly) 




(237) 




The nondimensional gains kjj are defined as in equations (210). 

The governing equation for the general solution for the position error yields the 
relation 

[s + itfj s 4 - A^kj x + k^y^s^ + + kj x k]y^s^ + X^k 2 y - ^ 2 x ^s - X^k 2 x k 2 yj 0 !(s) = 0 (238) 

Once again, the attitude equation contains a purely imaginary root, which leads to neutral 
stability of the position error. The trends of the spacecraft's motion for the rate plus 
rate integral control law are thus similar to those for the pure rate control law. 
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The stability regime for this control law is shown graphically in figures 27(a) 
and 27(c). The nondimensional gains k^ x , k^y, k 2 x , and k 2 y must now be selected 
to satisfy the conditions of equations (237). The first two of these conditions are identi- 
cal to those for the pure rate control law and are given by figure 27(a). The next two 
conditions lead to stable motion in the second quadrant of figure 27(c), subject to the 
last restriction which represents a compatibility relation between the rate and rate 
integral gains. This compatibility condition yields a hyperbola with the equation 


k 2 y 2 - 2 + ( 


k lx + k ly)' 


k 2x k 2y + k 2x + 


( k 2y ~ k 2x 


i 


klx + kly)^ + 


k lx k ly) 


= 0 


(239) 


as sketched in figure 27(c). The stable region in this figure is then the area between the 
upper segment of the hyperbola and the k 2 X sgn X and k 2 y sgn X axes. 

If the nondimensional rate gains are small in comparison with unity, as for most 
practical control systems, then equation (239) can be approximated by 


( k lx + k ly) ( k 2x " k 2y) 


(240) 


and the resultant stability characteristics are given in figure 27(d). The spacecraft's 
motion is stable if the rate integral gains are selected from a triangular area of the 
fourth quadrant for a set of stable rate gains. 

By examining equations (235) and (236), it can be noted that the modified charac- 
teristic equation for the rate error reduces to a cubic equation when either of the two 
rate integral gains vanishes. This special case would occur during single-axis control 
of the spacecraft and is thus of particular interest. 


Stability restrictions are given by the standard rate gain restrictions and the 
relations 


0 > k 2 x sgn X > 
0 > -k 2 y sgn X 


( k lx + k ly)( 1 + ^ly)** 11 O 



(241) 


The resultant gains will fall on the boundary of the stable region in figures 27(c) 
and 27(d). Accordingly, stable control of the spacecraft is possible with single-axis 
rate plus rate integral commands. 

Numerical integration of the governing equations (235) is necessary to determine 
the rate histories for particular values of the control gains. Some general conclusions 
can, however, be drawn when the nondimensional rate gains are chosen to be consider- 
ably larger than the nondimensional rate integral gains. Inspection of equation (236) 
shows that the last two terms in this equation are now small during an initial transient 
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period after a disturbance, when the rates are large and their integrals are small. The 
rate equation is thus approximately equal to equations (217) and the rate gains can be 
selected from equations (226) and (227) to yield the desired damping characteristics 
during the transient period. As steady- state conditions are approached, the rate integral 
terms will predominate and the resultant control torques will tend to eliminate any resid- 
ual rate errors. The net effect of this law is a reduction in the gain magnitudes since 
the high rate error gains are no longer needed to reduce standoff errors in the steady- 
state condition. Damping and attitude hold characteristics for the rate plus rate integral 
law should thus be quite efficient. 

Rate plus attitude control law: Next, consider control torques developed from both 
rate and attitude errors and given by 


g X i x ( K lx n x + ^3x^) K k lx fi x + ^Sx^) 
gy = J“( K ly s V + K 3y^) = *( k ly% + J 


(242) 


The resultant equations of motion are 




f2x - + l-j — jfty - aAkg x e - f x 




'Mxl 


Oy - ~j “ a ^3y^ — fy 


(243) 


and yield the general characteristic equation for the complex rate error 

A 


X^kix + k ly ^ s 3 + o 2 + X 2 (l + ki x kjy^ s 2 + a 


2 2 
k 3y ^x*x k 3x VV] 


I* 


- oX 


( k lx + k ly) 


+ a 2 X 2 jk lx k ly + k lx - kiyj = 0 (244) 


The associated Hurwitz stability criteria give the following conditions for stability of the 
rate errors: 
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-Xkjy > Akj x 


k lx k ly > 


1 + ' X 


ksy^A x I x kg x ^AyI^ 


> crA 


( k lx + k ly) 


^ k 3y^^x^x k 3x^^-ylyj > \/ k lx + k ly 


o 2 - A 2 (l + k lx k ly ) 


I “ a - 2 ( 1 + k lx k ly) + 4o 2 A 2 (l + kjy - ki x ) 


(245) 


After substitution of equations (215) into equations (242), the characteristic equa- 
tion associated with the general solution for the attitude error becomes 


- A|ki x + kiy^s^ + o 2 + X ^1 + ki x kiy| s 2 - c^A 


I V + I V " 1 9 


( k 3y k 3x 


L X Xy - X Z )[ — - - 


+ (k lx + k ly ) s + c^A 2 |l + (k lx - k 3x )(k ly + k 3y ) 


x *y 

= 0 (246) 


Corresponding stability conditions for the attitude errors are 


-Akiy > Aki x 


"\ 


r / \ 2 

k lx k ly > - 1 + 


>( > ; * !? V 


l y !x 


y J X + ly ” I 2 


( k lx " k 3x)( k ly + k 3y) 


> -1 


> (247) 


c _3y 

lx 


^ ( k lx + k ly) ( k 3x k 3y “ k lx k 3y " k ly k 3x) > (*x + ty “ ^z) 

- w) + V - ■ if) + (“lx * k ly) [«* - k2 (l ♦ k lx k ly^ 
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Both equations (245) and (247) must hold for stable error histories. Usually, the 
nondimensional rate and attitude gains are small in comparison with unity, and the sta- 
bility conditions given by these equations can be approximately represented by 


"X 

-Xk ly > Ak lx 


k lx k ly > - 1 


0 > 



k 3x \ 
lx + k lyy 



c 3y 


y klx + ^ly J 





> 1 


J 


(248) 


The rate gain restrictions are now identical to those for the pure rate law and are thus 
shown in figure 27(a). The attitude gain restrictions, illustrated in figure 27(e), yield a 
stable region falling between the lines 


k 3y 


k lx + k ly jVz) k lx + k ly/\Iz J 


k 3y 



‘3x 


"3x 


k lx + k ly/\ I z/ \ k lx + k ly/\ I z, 


2|-1| 


(249) 


It is worthwhile to note that either of the two attitude gains may be zero, and attitude and 
rate damping of the spacecraft is thus possible with torques applied about a single space- 
craft axis. The error signals needed in the calculation of these required torques can be 
taken from a single rate gyro and a single sun sensor mounted on that axis, and the 
mechanization of such a control system appears to be very simple. 

As seen from equations (244) and (246), the determination of the damped rate and 
attitude errors again requires the numerical integration of equations (243). It is 
observed, however, that both damping and reorientation control can be provided by the 
rate plus attitude control law and that this law can correspondingly be used to maneuver 
the spacecraft. In comparison, the rate and the rate plus rate integral control laws were 
restricted to holding an already established inertial position. 

In addition, it may be noted that pure attitude control, for which kj x = k^ y = 0, 
will result in several zero coefficients in equations (244) and (246). The associated 
spacecraft motion is, at best, neutrally damped and may diverge for certain forcing 
functions. Pure attitude control then provides no damping of the spacecraft rate and 
attitude errors. This conclusion is, of course, in agreement with previous results. 
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General considerations: Other control laws may be investigated in an identical 

manner by selecting the control torques, developing the complex governing equations, 
and defining the resulting stability regions for the spacecraft's motion. If the control 
system gains are chosen to satisfy these criteria, the motion of the spacecraft will be 
damped. The determination of time histories of the controlled spacecraft motion 
becomes a rather tedious task, however, and is perhaps easiest if the governing equa- 
tions (212) and (215) are programed on a digital or analog computer. After it has been 
established that the selected control gains lead to stable motion of the spacecraft, par- 
ticular disturbance functions can then be studied on the computer to determine the 
resultant spacecraft time histories. 

Extension of the method to include nonlinear control commands is possible, but 
nonlinear techniques (refs. 21 to 27) must then be used to define the stability of the 
governing equations. A preliminary selection of the type of on-off control commands, as 
represented by step torques, may be made by noting that the spacecraft response for an 
amplitude- limited control system with high gains approaches that for an on-off control 
system. Sensor inputs and the signs of the control torques may thus be chosen from the 
proportional analysis. Time histories for the on-off commands can then be obtained by 
substitution of the corresponding step functions in equations (42). The solutions for the 
body rates and Euler angles now are found by a piecemeal process, and the forcing and 
solutions functions change whenever the dead bands for the on-off system are crossed. 

Linear control optimization (ref. 28), as represented by a minimum mean square 
error criterion, may also be considered. Maximum torque or error limitations can be 
included as restraints in such an analysis. 

Control sy stem selection .- Having determined a control law which leads to accept- 
able damped motion of the spacecraft, physical systems that can develop the actuator- 
torque histories required by the control law must next be selected. The choice of such 
control systems is generally made on a minimum-launch-weight basis, and relations 
between the control system weight and its impulse or momentum storage capacity are 
needed to evaluate the comparative merits of various control hardware. Preferably, 
these relations should not necessitate the detailed design and optimization of competitive 
systems for a particular spacecraft. 

An empirical representation of the total control- system launch weight in terms of 
the angular momentum or impulse provided by the system is used in this analysis. Such 
a representation gives reasonable approximate values for the control- system launch 
weight and allows the rapid comparison of different control actuation schemes. Further- 
more, the empirical results are completely independent of the spacecraft inertia char- 
acteristics or dynamics. 
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Two tasks then remain, namely, sizing and implementation of the control system. 
To size the system, the spacecraft’s angular momentum envelope is determined by inte- 
gration of the torques corresponding to simultaneous application of all "worst-case" 
disturbances. The launch weight for the control system can then be developed from the 
empirical data, and a preferred system concept can be selected. The implementation 
of this concept requires the solution of the control system equations to define the actual 
torques the system must generate in order to provide the desired control law and elimi- 
nate cross coupling moments caused by any angular momentum stored in the control 
system. 

Control system components may be divided into the general classes of momentum 
storage units and reaction control units. Momentum storage units comprise reaction 
wheels, single-gimbal control moment gyroscopes, and double-gimbal control moment 
gyroscopes. Reaction control units consist of reaction jets with variable mass flow or 
pulse modulation. 


The angular momentum envelopes for these two classes of control components are 
then given by 


h M = h Mx + iH My 


for the momentum storage units and 


h R = H Rx + iH Ry 


~ ^ — 1 
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( f x at 
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0 M dt 

+ i 

n t<J 

%J 0 M dt 


(250) 


(251) 


for the reaction control units. The integra- 
tion is carried on over td, the time interval 
of application of the "worst-case" 
disturbances. 

Reaction wheels: For sizing purposes 

the reaction wheel is taken to be a flywheel 
which is accelerated by means of a torque 
motor to produce reaction torques on the 
spacecraft. A sketch of a control system 
with two such wheels is shown in figure 28. 



Figure 28.- Spacecraft control with reaction wheels. 
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From manufacturer's data, the basic weight of a reaction wheel (ref. 29) with a 
minimum alternating- current motor configuration is 


W w = 6.3 + 170 


H 


W 


Wt 


(252) 


where W W is the total reaction wheel and motor weight in lbm, Hy/ is the angular 
momentum capacity in ft-lbf-sec, and o>g is the synchronous wheel speed in rad/sec. 

Similarly, the reaction wheel power may be derived from empirical data (ref. 30) as 


P W ~ 2.77GgCOg 


(253) 


where Gg is the stall torque in ft-lbf. If a power weight conversion factor of a 
lbm/watt is introduced, then the equivalent power system weight for the reaction wheel is 


W p = 2.77aG s w s 

and the total weight chargeable to one reaction wheel becomes 


H 


W 


W wx = W w + W p = 6.3 + 170 —2. + 2.77aG s w s 


(254) 


(255) 


To optimize the total weight for a given angular momentum and stall torque, the total 
weight is differentiated with respect to wheel speed and the result is set equal to zero. 
Substitution of the corresponding wheel speed into equation (255) gives 


WwT = + 21.7yaGgH^y 


(256) 


The control torques may be assumed to be sinusoidal with amplitude Gg and fre 
quency X so that 

iXt 


g = -G s e 


(257) 


can be taken to be a good approximation to the control moments. Since these control 
moments are equal to the total rate of change of the angular momentum components for 
the reaction wheels, it follows that 

Gg cos Xt = - a H\v2l 


> 


or by introducing 


Gg sin Xt = H W 2 + <*H W1 


^xy ~ H W1 + ^W2 


(258) 


(259) 


^sed by courtesy of the Bendix Corporation. 
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that 


( 260 ) 


■j f ■ 

Gg e = H x y + iCHjfy 


Expression (260) can be integrated to give 


G c 


Hvv = 1 -J*L)e iXt 


(261) 


ix y ~\o + a, 

for no initial wheel momentum. Maximization of equation (261) further yields the value 


H _ G s 

for each of the two reaction wheels. 

Substitution of equation (262) into equation (256) leads to the expression 


( 262 ) 


W wx = 6.3 + 


21.7\/a(cr + A) H w 


(263) 


which is plotted in figure 29. 

The weight of the reaction wheel control system is now given by 


W T - W wx (H Mx ) + W wx (H My ^ 


(264) 


where W wt (H M x ) and 
W WX ( H My) are obtained from 
figure 29 with values of Hj^ x 
and HMy determined from 
equation (250). 

Note that the spin rate a 
must be very small if reaction 
wheels are to be efficient. As an 
example, a power conversion factor 
of 1 lbm/watt and a spin rate of 
0.25 rad/sec for a flat disk config- 
uration would yield W x = 3074 lbm 
for an angular momentum require- 
ment of 100 ft-lbf-sec along each 
axis. Since such exorbitant weight 
penalties are impractical, reaction 
wheels generally are inacceptable 
for the damping control of spinning 
spacecraft. 



Figure 29.- Variation of total reaction wheel weight with required 
angular momentum. (Conversion factors for SI Units are given 
on page 10.) 
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If such wheels are used for spacecraft with very low spin rates, the governing 
torque relations become 

g = -G = -H xy - icrH xy (265) 

and 


T z = OyH x - fi x H y (266) 

where G is the complex torque applied to the reaction wheels and T z is the cross 
coupling torque applied to the spacecraft by the control system. Since the body rates 
are small, this sinusoidal cross coupling torque is relatively small and its effect on the 
spin rate a can be neglected. 


The desired control torques are then obtained by directly applying the reaction 
wheel torques 

G = ^W1 + = "S (267) 

by means of the wheel actuators. 


Control moment gyroscopes: A control moment gyroscope consists of a flywheel 
which spins at a constant speed and is mounted on a single- or double- gimbal arrange- 
ment. Control torques are now developed by precessing the flywheel. Torque actuators 
mounted on the gimbals provide the necessary precession torques. Sketches of control 
systems using single- and double-gimbal gyros are given in figures 30(a) and 30(b), 
respectively. 



Figure 30.- Spacecraft control with control moment gyros. 




Since the weights of single- and double-gimbal gyros do not differ appreciably, 
both these units are assumed to have identical launch weights. The basic weight of a con- 
trol moment gyroscope can again be developed from manufacturer's data (ref. 31) and 
becomes 

W G = 1.37H g 0 - 68 (268) 

Power requirements now are derived primarily from the windage and friction losses for 
the flywheel and can be approximated by empirical data derived from computer analyses 
(ref. 32) as 

Wp = 1.47aH G °* 362 (269) 

The power required by the gimbal actuators is small and is neglected. The total weight 
of the gyro is then 

Wqt = W G + Wp = 1.37H g 0 - 68 + 1.47aH G 0 ' 362 (270) 

Launch weights of the gyro are plotted against angular momentum in figure 31, and a 
comparison with figure 29 shows that the total weight for a system using gyros is much 
less than that for a system using reaction wheels. 



Figure 31.- Variation of total gyro weight with required angular 
momentum. (Conversion factors for SI Units are given on 
page 10.) 
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In a comparison of single- and double-gimbal control moment gyros, it is observed 
that the weight for the single-gimbal gyro system is derived from 

W T = w Gt( h Mx) + w GT( H My) ( 271 ) 

whereas the weight for the double-gimbal gyro system becomes 


w T = w GT ! 


H- 


Mx 


+ H 


My 


cos 


Og 


(272) 


where ctg is the limiting gimbal angle and Hj^ x and Hj^y are again found from 
equation (250). Since o tg is generally 60°, the double-gimbal gyro system is somewhat 
lighter than a control system using two single-gimbal gyros. 

To derive the gimbal actuator commands for the implementation of the desired con- 
trol laws, these torques are again deduced from the total rate of change of the angular 
momentum vector for the gyro system. Thus, the complex control torque is 


g — Gxy — Hxy - i joHjjy - H z S2x^J 
and the cross coupling moment becomes 


(273) 


T z = -H z + OyHjj - OxHy 


(274) 


The minor changes in the spin rate a due to T z are neglected in this linear formula- 
tion, and the necessary control commands are now found by expressing Hxy and Hz 
in terms of the gimbal angles. 


For the single-gimbal gyros this gives 

H 


X y — Hx + iHy — Hqi sin 0g — sin 


H Z — Hqj cos 0g + Hq2 cos •Pg 


(275) 


The governing torque equations now reduce to 


g - -Gjjy - - 


H G -j^0g + fiy^cos 0g + Hg2( ct sin <pg + f2y cos 


+ l 


and 


i H G2 ^g - S^x^cos cpg + H G i^c sin 6g - Ox cos 9gj 


T z — H G j^0g + Oj^sin 0g+ H G 2^g - Ox^sin cpg 


(276) 


(277) 
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where the small gimbal accelerations have been neglected. Gimbal torques may be 
commanded directly from 

Gxy = G x + iGy = -g (278) 


in an open- loop system. 

For the double-gimbal gyro the angular momentum components along the space- 
craft's axes become 


Hxy = HQ^sin 0g - i sin cpg cos 6gj 


and 


H z = Hq cos (fig cos 6g 
The torque relations thus are 

g = -Gxy = -^G x + iGy cos (Pgj 


= -HG|j^g + cr sin (pg + S2y cos <pg cos 6g - i + STxjcos 6g cos 
- (o + dg sin cpgjsin 0gjj> 


cpg 


and 


where the gimbal acceleration terms are again neglected. 
Necessary gimbal torques are developed from 

Gx - -IxSx 

Gy = -Iygy sec cp g] 


(279) 


(280) 


T z = -Gy sin cpg = H G |^(pg + £2x)sin cpg cos B g + (py + 6g cos <p g )sin (281) 


(282) 


and require the measurement of the gimbal angle cp g . 

Reaction jets: The reaction jet system comprises the propellant, oxidizer, engines, 

and tankage weight necessary for the spacecraft's control. To arrive at weight estimates 
for such a system it was assumed that the usable specific impulse considering engine effi- 
ciency, expulsion efficiency, and ullage would be 290 lbf-sec/lbm and that the propellants 
would be storable hypergolics housed in tanks with positive feed expulsion diaphragms. 
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Manufacturer's data^ can then be extrapolated (ref. 33) to yield the idealized total system 
weight 

W JT = 0.01011-p 0,912 (283) 

which is presented in figure 32. Here I T is the total impulse in lbf-sec. This impulse 
may be written in terms of the total momentum envelope for the spacecraft as 


I T = 


H 


Rx H Ry Hrz i 

Zy Zy If 


where 


ly y and 


£. X ty /r£ 

are the moment arms in ft about the 


(284) 

X, Y, and Z body axes. 



Figure 32.- Variation of total reaction jet system weight with required total 
impulse. (Conversion factors for SI Units are given on page 10.) 


2 


Used by courtesy of the Minneapolis- Honeywell Regulator Company. 
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An assessment of the weight of the reaction jet system thus requires the develop- 
ment of total momentum envelopes for the spacecraft mission. The momentum for par- 
ticular disturbances must be determined from equation (251) and the resulting momentum 
components along the spacecraft’s axes must then be multiplied by the probable frequency 
of occurrence of each disturbance. By repeating this process for all disturbances and 
summing the individual momenta along each axis, a total momentum envelope per 
sampling period is obtained. The weight crossover time between momentum storage and 
reaction control system can now be established from figure 32, since the total momentum 
envelope per unit time has been developed. If the mission time exceeds this crossover 
time, momentum storage systems should be selected for damping of the spacecraft's 
motion. 

If reaction-jet systems are chosen, they can be combined with a mass-balancing 
system which compensates for any constant products of inertia resulting from crew 
motion or cargo transfer and eliminates limit cycling of the jets about the new principal 
moments of inertia. Such a system could, for example, pump the propellant to different 
positions within the spacecraft to obtain its control torques. Since the design of this 
mass-balancing system is very much dependent on the spacecraft geometry, it is not 
considered here. 

Mostly, however, the control system consists of both reaction jets and momentum 
storage systems. The jets then provide for attitude control and orbit keeping and the 
momentum storage system is used to damp any oscillatory motion of the spacecraft. 
Attitude control commands are now used to actuate the reaction jets while rate and rate 
integral commands provide control laws for the momentum storage system. The develop- 
ment of such combined systems is again dependent on specific spacecraft and disturbances 
and will not be attempted in this analysis. 

Comparison of Exact and Approximate Solutions 

Two possible manned spacecraft were considered for a comparison of the results 
of the numerical integration of the exact equations of motion and the results of the pres- 
ent analytical solution. These spacecraft were a cylindrical manned orbital research 
laboratory and a large hexagonal space station. 

Manned orbital research laboratory .- The manned orbital research laboratory 
(MORL) is proposed as an earth-orbital laboratory in which scientific and engineering 
experiments could be conducted over extended time periods. The basic laboratory is 
designed to support a crew of six astronauts in a 200-nautical-mile orbit for periods up 
to 5 years. During spinning operation, the laboratory module and the last stage of its 
Saturn booster would remain attached by a system of cables and would rotate about a 
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common mass center. The resultant centrifugal force would produce an effective gravity 
field in the manned module. A sketch of the corresponding MORL configuration is shown 
in figure 33. Assumed inertia and mass characteristics (ref. 34) for this configuration 
are given in table VI. The inertia distribution is near- cylindrical, and the manned 
module and booster counterweight rotate about the Z-axis at 0.4 rad/sec. A gravita- 
tional acceleration, equal to one-fourth that at the earth's surface, acts on the manned 
module because of this rotational rate. 

The effects of various disturbances on MORL are summarized in table VII. The 
disturbances include residual rate and attitude errors after spin-up, moments applied by 
an attitude control system valve failure in the open position, and several "worst-case" 
crew motions. These motions comprise step translation to an extreme position within 
the laboratory and linear oscillations which could result from trampoline exercise, lad- 
der climbing, or floor pacing. A linear velocity of 4 ft/sec is selected for all oscillatory 
motions, and the entire crew of six is taken to be a single equivalent mass with a mass 
factor Q of 36 slugs. 

Equations defining these disturbances are listed in the second column of table VII. 
The resulting error limits have been found by hand calculations of the analytical upper 
bounds and by extrapolation of the error time histories obtained from numerical integra- 
tion of the exact equations of motion on an IBM 7094 system. 



L-67-966 

Figure 33.- Artist's sketch of possible manned orbital research laboratory. 
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Both rate and attitude error limits are given in the table. The rate error limits 
range from about 0.004 rad/sec for the step inertia products to about 0.024 rad/sec for 
the step torques, and the attitude error limits vary from approximately 0.01 rad/sec for 
the step inertia product to approximately 0.18 rad/sec for the residual errors after 
spin-up. Significantly, the errors caused by the periodic mass motions are several 
times greater than those produced by instantaneous motion to a final position. Predic- 
tions of maximum spacecraft errors due to crew motion must thus consider any periodic 
crew motions that may occur. Onboard experiments, which require high-accuracy con- 
trol of the spacecraft, may be adversely affected by the oscillatory crew motions and 
may require restriction of these motions. 

The approximate error limits developed from the analytical solutions show reason- 
able agreement with the maximum errors determined from the exact solutions. The 
deviations of the approximate error limits from the exact error limits are generally less 
than 20 percent of the exact error limits. These upper bounds of the spacecraft errors 
thus give a conservative estimate of the effects of various disturbances and should be suf- 
ficient for initial engineering design applications. 

Two of the solutions described in table VII have been selected for a comparison of 
the actual error histories and the error histories given by the analytical solution. The 
disturbances are the step product of inertia and the vertical mass oscillation. Both the 
uncontrolled and the controlled spacecraft motions were considered. Solutions for the 
uncontrolled motions were obtained from numerical integration of the exact equations of 
motion and from evaluation of the error relations developed in this analysis. Solutions 
for the controlled motions were obtained from numerical integration of both the linearized 
and the exact equations of motion. All calculations were carried out on an IBM 7094 
system. 

Single-axis control commands, which apply torques about the X or minimum- 
inertia axis of MORL, are most efficient. The nondimensional rate damping gain kjy 
in equation (227) may accordingly be set equal to zero. This operation yields 

k lx = - = - 2r D ( 285 ) 

as the nondimensional rate damping gain for the X-axis. The corresponding rate integral 
gain, when used, is arbitrarily selected as 



to fall in the stable region of figure 27(d). 
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For MORL, the time constant t D is equated to one spin cycle or 20 sec. Values 
of the control constants are then 

Kq x = -13 300 ft-lbf-sec/rad 

and 

K 2x = -2 040 ft-lbf/rad 

The rate damping ratio becomes 

r D = 0.1 

and state-of-the-art gyros (ref. 35) in the 500 to 1000 ft-lbf-sec class (fig. 31) and jet 
hardware can provide the necessary control torques. 

The MORL response to the step product of inertia is illustrated in figures 34 
and 35. The uncontrolled rate and attitude errors are given in figure 34. It is apparent 
that the exact and approximate solutions are virtually identical. Both the rate and the 
attitude errors are biased sinusoids. As expected from the trace analysis, the largest 
rate error occurs about the Y-axis and the largest attitude error corresponds to rotation 
about the X-axis (figs. 9 and 10). To the crew, the rate error appears as a minute 
rolling motion of the laboratory floor with a maximum amplitude of 0.6°. 





e, 

deg 



L 

0 


20 


_! L 

40 60 

t, sec 


-I 

100 


l _J ! ! I I 

0 20 40 60 80 100 

t, sec 


Figure 34.- Uncontrolled laboratory error histories for step inertia products. 
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Figure 35.- Laboratory error histories for step inertia products and pure rate control. 
(Conversion factors for SI Units are given on page 10.) 


The controlled response of the laboratory to this disturbance is shown in figure 35. 
This figure, which corresponds to pure rate control about the minimum inertia axis, 
again gives the same results for both the exact and the approximate solution. The labo- 
ratory oscillation is reduced to steady coning in about three spin cycles. The constant 
residual rate errors produce constant control torques which counteract the mass unbal- 
ance torques produced by the products of inertia. The effective disturbance torque is 
thus less than it was for the uncontrolled response, and the residual rate errors are cor- 
respondingly somewhat smaller than the constant components of the uncontrolled rate 
errors. The oscillatory terms in the uncontrolled rate errors are due to the accelera- 
tion terms associated with the introduction of the inertia products and tend to zero in the 
controlled error histories. 


The attitude errors for the damped rate errors become 

-iat ^xyr 
Q! r - <p v + i0 r - o^e — — 


(286) 


from equations (215). For the step products of inertia, the attitude errors result from 
changes in the body rates which do not affect the total spacecraft momentum. The con- 
tribution Oq of the transient oscillatory terms in the body rates to the attitude errors 
tends to zero, and residual attitude errors are given by 



(287) 
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where S2x r and Oy T denote the residual body rates. Both the rate and attitude errors 
approach constant values for the theoretical solution. 

As predicted by the analysis of the controlled spacecraft characteristics, single- 
axis rate control is acceptable for normal operation and experiments which do not require 
high-accuracy stabilization of the spacecraft. The MORL response to a vertical periodic 
motion of the entire crew is depicted in figures 36 and 37. Figure 36 illustrates the 
uncontrolled results. The exact and approximate solutions check very closely. The rate 
and attitude errors now comprise a low-frequency, large- amplitude sinusoidal oscillation 
due to precession within the outer (X) ellipses and high-frequency, small- amplitude oscil- 
lations due to precession within the inner (p) ellipses. (See figs. 12 and 13.) The maxi- 
mum errors are two to three times as high as the corresponding errors for the step 
inertia product. The laboratory floor also undergoes irregular rolling motions with maxi- 
mum amplitudes of about 2°. Since the distance (ref. 34) from the center of rotation to 
the laboratory floor is approximately 50 ft, this roll can produce a 2-ft total translation 
of the station floor and could present some difficulties to a moving astronaut within the 
laboratory. 
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Figure 36.- Uncontrolled laboratory error histories for vertical mass oscillation. 
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Figure 37.- Laboratory error histories for vertical mass oscillation and rate plus rate integral control. (Conversion factors 

for SI Units are given on page 10.) 


The controlled laboratory motion is presented in figure 37. Control torques are 
applied about a single axis, but a combined rate and rate integral control law is used. 

The laboratory motion is quite similar to the uncontrolled motion but exhibits damping of 
the free vibrations, as is apparent from the gradual decrease in the corresponding error 
terms. Since the periodic forcing function terms predominate, this type of response is 
to be expected. Agreement between the approximate and exact solutions is very good, 
and the small differences in the error histories can only be detected for the angle cp. 

The uncontrolled and controlled error histories developed from the approximate 
solution for the nonsymmetric MORL are practically coincident with the exact error 
histories for all the disturbances that have been examined. The analytical solution is 
thus a useful tool for the study of the nonsymmetric laboratory motion. 

Large manned space s tation .- A second possible manned rotating spacecraft is the 
large spinning space station, such as the 150-ft station which is considered here. This 
station (ref. 1), shown in figure 38, has six cylindrical outer modules arranged in the 
shape of a hexagon. The outer modules are connected to a central hub and docking port 
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Figure 38.- Artist's sketch of possible 150-ft manned space station. L-67-967 


by three spokes. Rotation about the maximum inertia axis provides artificial gravity for 
the living modules. The crew of this space station would vary from 6 to 21 astronauts. 

Assumed characteristics for the 150-ft space station are listed in table VE3. The 
inertia distribution approaches that of a flat disk and the spin rate is 0.314 rad/sec. The 
crew is assumed to be composed of six astronauts with an effective mass factor Q of 
about 36 slugs. 

Disturbance effects on this space station are summarized in table IX. The dis- 
turbances are similar to the MORL disturbances, and periodic crew motions are simu- 
lated by motion of a single equivalent mass with a linear velocity of 4 ft/ sec. The rate 
error limits range from approximately 0.006 rad/ sec for the step inertia products to 
approximately 0.015 rad/sec for the circumferential mass motion. The attitude error 
limits vary from about 0.017 rad/sec for the step inertia products to about 0.18 rad/sec 
for the residual errors. The errors due to periodic motions are considerably greater 
than those introduced by the step products of inertia. 

Circumferential crew motions and the residual errors were chosen for a further 
comparison of the approximate and exact solutions. Both uncontrolled and controlled 
solutions were developed. 


105 



Control torques were now applied about both station axes and the corresponding 
control gains were assumed to be equal. Equations (226) and (227) now yield 

k lx = k ly = “ = " 

as the nondimensional damping gains. The corresponding rate integral and attitude gains, 
when used, are selected as 

k lx 

k 2x = _k 2y = — 

and 

k 3x = -k 3y = k lx 

from the stable regions of figures 27(d) and 27(e). The physical control gains become 

Klx = Kfy = -222 817 ft- lbf- sec/rad 

K2 x = -K 2y = -75 000 ft- lbf/ rad 

and 

k 3x = -K 3y = -70 000 ft- lbf /rad 



from equations (210). A time constant of about three spin cycles or 54 sec was selected 
to give the damping ratio 

r^ = 0.02 

Higher values of damping would require exorbitant control moment gyro and reaction jet 
control systems. Even the selected value requires gyros in the 5 000 to 10 000 ft-lbf-sec 
class and exceeds the present state of the art in gyro hardware shown in figure 31. Rapid 
jet damping, although feasible, results in large fuel consumption. 

The advantages of single-axis control for nonsymmetric vehicles become obvious 
when it is noted that the MORL, with about one-half the spin momentum of the 150-ft 
station, requires a control system that is an order of magnitude smaller. In addition, 
the MORL is able to achieve lower time constants and considerably better damping ratios. 
These results lead to the conclusion that nonsymmetric spacecraft, spinning about a maxi- 
mum inertia axis, are preferable from the control standpoint and that single-axis stabili- 
zation about the minimum inertia axis can result in a major saving in control system 
weight for these spacecraft. 
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The 150-ft- station motion for the circumferential mass transfer is given in fig- 
ure 39. The approximate and exact solutions are in good agreement, and the time 
histories exhibit slow oscillations. These oscillations (see figs. 21 and 22) consist of a 
large-amplitude sinusoid with the mass motion frequency p and small- amplitude sinus- 
oids with, approximately, the precession frequency X. The angular deviation of the 
gravity vector has an amplitude of about 3° and appears as a corresponding slow rolling 
motion to the crew. For the 150-ft space station this rolling motion produces a 4-ft 
oscillation of the station floor. 
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Figure 39.- Uncontrolled space station error histories for circumferential mass motion. 


The controlled space station response is illustrated in figure 40. Control is 
derived from rate plus rate integral commands, and the station rates are damped to a 
purely sinusoidal trace in approximately nine spin cycles. The constant rate term in the 
damped trace eventually disappears under the action of the rate integral commands. 
Residual rate errors may then be expressed as 

ftxy r = -f2 r e ipt (289) 
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Figure 40.- Space station error histories for circumferential mass motion and rate plus rate integral control. (Conversion 

factors for SI Units are given on page 10.) 

where fir denotes the half-amplitude of the residual rate. From equation (160), this 
residual rate amplitude term is 

„ „W a + p ) 2 

W - P) 


so that 


^xyi 


I rz (a+ p)‘ 


ipt 
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The corresponding residual rate error becomes 
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I r z( CT + P) 

V a + P/~ 

Li<*- p) J 


cos pt 


The residual terms correspond to the coefficients of e*^ in the uncontrolled solution 
functions. The control system thus has little effect on the magnitude of the errors 
directly caused by the constant circumferential mass motion, as would be expected for 
the selected low value of the damping ratio. The approximate solution compares favor- 
ably with the exact solution for this example. 


Figure 41 presents the space station response for residual rate and attitude errors. 
As anticipated from figure 15, the rate errors are simple sine and cosine curves. The 
attitude errors, defined by equation (133), are somewhat more complex sinusoids. The 
exact and approximate solutions checked to within three significant figures. 
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Figure 41.- Uncontrolled space station error histories for residual errors. 

The controlled space station motion with the residual errors is shown in figure 42. 
Rate plus attitude control commands are now employed, and the station completes the 
required 10° reorientation about two axes in approximately 16 spin cycles. The analyt- 
ical and exact solutions again were identical. 
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Figure 42.- Space station error histories for residual errors and rate plus attitude control. (Conversion factors for 

SI Units are given dn page 10.) 

Summary of comparison .- The analytical solution was in excellent agreement with 
the exact solution for all examples considered. Since the spacecraft used in the com- 
parison are typical of future rotating manned spacecraft, the analytical solution should 
be valid for the determination of the dynamics and control of most such spacecraft. Ana- 
lytical results for unmanned spacecraft, which may have larger torque disturbances and 
residual errors but have few or no inertia changes, should also be acceptable. The time 
history data obtained for the step torques and residual errors were accurate to three 
places and an increase in these disturbances should not appreciably degrade the results 
in the linear range. The analytical solution thus offers a simpler, more economical, and 
more direct means of assessing the effects of various disturbances and spacecraft char- 
acteristics on the spacecraft motion than the computer runs. The insight into the 
mechanics of motion, which is gained from the error formulation developed in this analy- 
sis, should be of major value to future work on the dynamics of arbitrary rotating 
spacecraft. 
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ANALYSIS OF SPIN-UP AND DESPIN MODE 


Governing Equations 

For the present application of the governing equations for the spin-up and despin 
mode, the spacecraft disturbances are assumed to be restricted so that no internal mass 
movements occur and no moments are exerted about the spacecraft X- and Y-axes. 

With this restriction, equations (8) and (12) then become 


and 


Q, z 


Iz 


^zO^zo 



(290) 


^ = ^O^zo 



M z dt 


dt 


(291) 


Spin-up and despin moments about the Z-axis are assumed to be provided by 
constant-thrust, pulse-modulated jets (refs. 1, 34, and 36). Since the control of the space- 
craft during this mode is quite straightforward, the main problem is the selection of a 
spin-up and despin technique which minimizes the associated fuel consumption for rigid 
and extensible spacecraft. 


Rigid Spacecraft 

For the rigid spacecraft configurations, such as the large hexagonal space station, 
the jet moment arms remain constant. Spin-up and despin fuel is thus given by 


w = ^z*f _ *zffozf 
SU ^zf I SP fzf I SP 


(292) 


where the subscript f denotes conditions after completion of the spin-up maneuver. 
The simplest spin-up technique would apply continuous thrust or constant -width thrust 
pulses until the desired spin rate D z f is reached. 


Extensible Spacecraft 

Math ematical model .- For extensible spacecraft configurations, such as the MORL, 
the fuel calculation becomes somewhat more difficult. As an example, consider the 
mathematical model in figure 43. The spacecraft consists of a manned module with 
mass m m and a counterweight module with mass m c . The two modules are connected 
by a flexible cable or strut arrangement, which is extended to produce a large rotational 
radius. The distance between the module mass centers is designated l, the offset dis- 
tance between the thrust Pm and the manned module mass center is Z m , and the offset 
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Figure 43.- Mathematical model for spin-up fuel calculations. 

distance between the thrust P c and the counterweight module mass center is l c . 

These offset distances yield a larger moment arm and may be required to maintain the 
spin-up thrust line normal to the line connecting the mass centers. To minimize oscil- 
lations of the manned module about its Z-axis, some type of rate damping should also be 
provided. Rate damping moments can be supplied by a small reaction wheel or passive 
dampers. The individual module oscillations about their respective mass centers are, 
however, neglected for the spin-up fuel calculations. 

Spin-up thrusts may be produced by jets on the manned module, by jets on the 
counterweight module, or by jets on both modules. The first method is preferable when 
the counterweight module mass exceeds the manned module mass, and the second method 
is preferable when the manned module mass exceeds the counterweight module mass. 

The third method may be used if jets are mounted on both the manned and counterweight 
modules; a pure couple about the spacecraft mass center can now be produced by 
selecting 

( 293 ) 



Spin-up thrusts would be simultaneously applied to both modules for this method. 

All three thrusting methods lead to effective moment and moment arm relations of 
the form 


l z = DjZ + D2 
M z = D3 ^z 


( 294 ) 
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where 


-f-V _ 111 lx Vs 111 

1 ~ (Pm + Pc)( m m + m c ) 

Pm^m + P C Zc > 

° 2= Pm + Pc 


D3 - Pm + Pe 


for the spacecraft model. 

The spin inertia I z can be written as 




- Iq _ + D gZ z 
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where 


Iq “ I m + Ic + 


' m c m m \ 
m m + m c y 



D 4 = 
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/ 2m c m m \/ d 2\ 

y m m + m cA Dl 2 J 

/ m c m m \ / ! x 

\m m + m cj( Dl 2) 




(295) 


(296) 


(297) 


and I m and I c denote the manned and counterweight module inertias, referred to the 
module mass centers. 

Equations (293) to (297) now allow the simple formulation of the total fuel consump- 
tions for different spin-up techniques. 

Spin-up an d despin techniques .- Spin-up is assumed to occur in the following man- 
ner. While rigidly coupled, the two modules are brought to an angular rate The 

modules are then separated by extending the flexible module connector under action of 
the centrifugal force. During this extension process the spacecraft momentum, spin- rate, 
or spin-up thrust may be held constant. After the full cable extension is reached, the 
two modules are spun up to the final spacecraft spin speed. Despin requires this 
sequence in reverse order. 
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The fuel required for spin-up or despin can be expressed by 



(298) 


where the subscript i denotes conditions after the initial spin-up to f2 z i, the subscript 
e denotes conditions after the extension, and the subscript f again denotes final 
conditions. 

Three characteristic spin-up techniques are considered here. These involve 
extension with constant momentum I z i^ z i, constant spin rate £2 z i, and continuous 
thrust. 

Constant- momentum extension: Perhaps the simplest spin-up technique is one 

where no spin-up thrust is applied during the extension. The modules are allowed to 
separate as desired while the spin speed automatically decreases to maintain the angular 
momentum constant. After achievement of the desired extension the spin-up jets are 
again actuated. 

The fuel consumption now becomes 


W SU = 


x zf fi zf 


Zz^SP 


1 + 


*zi^zi\ Azf 
*zf^zf/ Vzi 


(299) 


from equation (298). The rate of cgble (or strut) extension does not affect this fuel con- 
sumption and may be varied arbitrarily to maintain the cables in tension during the exten- 
sion. The fuel consumption is minimized by selecting the smallest value of O z i which 
yields sufficient cable tension at the completion of the extension. 

Constant- rate extension: For this spin-up technique, the spacecraft spin speed is 

maintained at its initial value f2 z i throughout the extension. The modules are again 
allowed to separate until the final extension is reached, and the spacecraft is then brought 
to its final spin speed f2 z f. 
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The fuel consumptions relation (298) reduces to 


W SU = 


^f^zf 

Zz^SP 



(300) 


for this technique. Fuel consumption now is minimized by selecting the smallest value 
of which yields sufficient cable tension at the beginning of the extension. 


The spin-up thrust for the constant- rate extension is established by the require- 
ment that the rate of change of angular velocity due to operation of the thruster must be 
greater than that due to the rate of extension or retraction. When the extension rate is 
maintained at a constant value Z z , this condition can be expressed as 


p 3 = ^zi^z 


2D P 


(301) 


The minimum thrust is thus directly dependent on the product of the extension rate and 
the initial spin rate. Corresponding conditions for variable extension rates may be 
developed from equation (290), if the time variation of l z is known. 

Continuous-thrust extension: Another possible spin-up technique would involve 
continuous thrusting during the spin-up. This brute-force technique requires rapid 
extension of the cable modules to be efficient, but is simpler to implement than the 
constant- rate extension. 

The associated fuel consumption is given by 


W SU = 


p 3tf 

X SP 



(302) 


To evaluate equation (302), a time history for either f2 z or l z during the extension 
must be selected. Values of O ze and t e can then be developed from this time history 
and equation (290), if it is recalled that l z - l z f and f2 z = J2 ze when t = t e . 


If the extension rate l z is constant, then 


fi ze ~ ^i^zi + p 3Zzi( t e " *i) + 2 (* e ~ *0 


“ *i - 


Zzf ~ Z ; 

lz 


Zl 


(303) 
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Substitution of equations (303) into equation (302) now yields 


TT7 _ ^Zf l . [ l 7l *Zi\ 


^zi^zi ^3(^zf ~ ^zi) 


f-zi 


22 , 


(304) 


and fuel consumption is optimized by selecting the smallest value of f2 z i and the largest 
value of i z which avoids cable slacking during the extension. 

Comparison of extension techniques: To compare the different extension tech- 
niques, the fuel consumption for each technique is expressed as the ideal fuel consump- 
tion at full extension plus an incremental fuel consumption for the extension process. 
From the ratio of these incremental fuel consumptions taken from equations (299) and 
(300), it follows that 


(^zi)cm t 

/ hi ) 

(° zi )cR,CT 

Wzi 2 / 



(305) 


is the condition corresponding to equality of the constant- rate and constant- momentum 
fuel consumption. Similarly the relation 


(^zi)cm 1 

( Iz i \_I 

°3 

fe. !) 

(^zi)cR,CT 

Wzi 2 / 2 

D 5^z( fi zi)c R> cT 

\ z zi / 


(306) 


must hold for equality of the constant-thrust and constant-momentum fuel consumption. 
The initial spin rates (fi z i) CT and (^zi)cR for the constant- thrust and 


constant- rate extensions should produce equal centrifugal forces to start the extension 
and are both assumed to be equal to the value (^zi)cR ct* s P* n ra ^ e (^zi)cM * or 

the constant- momentum extension is greater than or equal to (^zi)cR,CT* 

Equations (305) and (306) are represented graphically in figure 44. This figure 
allows the direct selection of the most economical spin-up technique for a particular 
spacecraft as a function of an extension length ratio Z z i/Z z f, a moment arm ratio T>2/lzi> 


a thrust ratio 


D 5^z(^zi)cR, CTj 


and a momentum ratio 


(^zi) 


CM 


_P zi )cR,CT J\D 5 Z z i 


A zi 
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Figure 44.- Fuel criteria for extension technique selection. 


The extension parameters D 2 , D 3 , l z , (ft z i) CR cx , and (^zi) CM are first chosen. 
The selected parameters and the initial spacecraft characteristics determine values for 
the moment arm ratio, the thrust ratio, and the momentum ratio. 

For a particular value of the extension length ratio, the moment arm ratio locates 
a fuel criteria point corresponding to the constant-rate extension and the thrust ratio 
locates a fuel criteria point corresponding to the continuous- thrust extension. The 
momentum ratio gives a third fuel criteria point for the constant- momentum extension 
which can be located on the ordinate. The lowest of the three fuel criteria values indi- 
cates the technique which yields the lowest fuel consumption. 


The special case, where the minimum centrifugal force for all three techniques is 
equal, is of interest since the minimum cable tension usually determines the extension 


parameters. For this case, 



( 307 ) 
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for the constant-momentum extension and 


D 3 = 




Da 

hi 

£2 

hi 



2 ^z(^zi)cR,CT 

Zzf 2 -^zi 2 


(308) 


for the continuous-thrust extension. 

From equations (307) and (308), 

Constant-momentum fuel criteria _ j + hi 
Continuous-thrust fuel criteria hi 


(309) 


and the continuous- thrust extension is now always more economical than the constant- 
momentum extension. Only the constant- rate extension and the continuous- thrust exten- 
sion need thus be compared. 

As an example, consider the manned orbital research laboratory described in 
table VI. The assumed extension parameters for this spacecraft are 

Zi = 37.5 ft 

Zj = 137.5 ft 

Mcr,ct ■ oa rad/sec 

£2 z f = 0.4 rad/sec 

D 3 = P m = 100 lbf 

D 2 = lm = 4 « 


and a spin-up technique which yields the lowest fuel consumption is desired. The cable 
tension must be greater than or equal to its initial value during the extension, and the 
spin-up time is immaterial. 

Since the minimum centrifugal force must be equal for all three techniques, the 
constant-momentum extension may be disregarded. From the given extension param- 
eters and equations (301) and (308), 

(^z)ct = 0.252 ft/sec 

and 

(z z ) CR = 0.172 ft/sec 

are obtained as the respective limiting extension rates for the continuous-thrust and 
constant-rate spin-ups. 
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The extension ratios now are 


and 

_P 5^z )ct ( nzi )cR, CT 

By using these values with figure 44, it is found that 

(Fuel criteria)^ = 1.0 

(Fuel criteria)^ = 1.6 

Since the fuel criteria for the continuous-thrust extension are considerably lower than 
those for the constant- rate extension, it follows that the continuous-thrust spin-up tech- 
nique requires the least spin-up fuel for this example. 

The actual fuel consumption values for the example, as computed from equa- 
tions (299), (300), and (304) for a specific impulse of 290 lbf-sec/lbm, are 

W SU,CT = 245 lbm 
W SU,CR = 256 lbm 
W SU, CM = 276 lbm 

The continuous-thrust spin-up requires approximately 12 percent more fuel than the 
ideal spin- up at full extension and the constant-momentum and constant- rate spin-ups 
require approximately 26 percent and 17 percent more fuel than the ideal value of 
219 lbm. A saving of about 22 lbm of fuel can thus be realized for each spin-up and 
despin cycle by selection of an optimum spin-up technique for this example. 

If desired, the response of the spacecraft to internal mass movements and external 
disturbance torques may be included in the analysis of the spin-up and despin mode by 
using equations (8) and (12) as the governing equations of motion. Supplementary linear- 
ized equations of motion for the relative module oscillations can be incorporated in such 
an analysis. 


tzi 
l zf 


= 0.334 


Do 

= 0.254 

Izi 
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CONCLUSIONS 


An approximate solution of the equations of motion of arbitrary rotating spacecraft 
with variable disturbance functions has been developed on the basis of small changes in 
the spacecraft body rates, Euler angles, and inertia terms. Complex representations 
have been used to define spacecraft and rate errors induced by the disturbance functions, 
and the solutions for the time history components and total error vectors have been 
examined for both uncontrolled and controlled spacecraft. 

The results of this analysis have led to the following conclusions: 

1. A comparison of the present analytical solution and solutions obtained by numeri- 
cal integration of the exact equations of motion for two typical manned spacecraft has 
shown that the analytical solution is in excellent agreement with the exact solution for 
the small angle and rate regime. The analytical solution provides a simpler, more 
economical, and more direct method of assessing the effects of various disturbances and 
spacecraft characteristics on the spacecraft motion and allows an insight into the mechan- 
ics of motion which cannot be derived from the numerical solution. 

2. Analytical upper limits of the rate and attitude errors induced by various dis- 
turbances are in reasonable agreement with the maximum errors found by interpolation 
of the numerical data. These upper limits should suffice for first estimates of the effect 
of the disturbances on the spacecraft motion. 

3. The spacecraft inertia distribution was found to have a significant effect on the 
spacecraft motion for equal disturbance characteristics. Spacecraft whose inertia dis- 
tribution approached that of a flat disk exhibit considerably more inherent stability than 
slender, near- cylindrical spacecraft spinning about a maximum or minimum axis of 
inertia. However, in practice the disturbance characteristics are directly related to the 
inertia distribution, so that the error bounds for spacecraft with different inertia dis- 
tributions tend to be similar. 

4. Periodic mass motions within the spacecraft may result in rate and attitude 
errors which are several times greater than those predicted for worst-case step prod- 
ucts of inertia. For equal disturbance characteristics, the largest errors resulted from 
circumferential mass motion in the direction of spin. Somewhat smaller errors were 
produced by radial mass oscillations in an offset spin plane and by vertical mass oscilla- 
tions parallel to the spin axis. Motions of the crew such as trampoline exercise, ladder 
climbing, or periodic translations along the spacecraft floor should be carefully examined 
to determine their impact on the spacecraft motion. 
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5. The spacecraft errors indicated instability trends when the spin axis became 
an intermediate axis of inertia during a mass motion and when the periodic motions took 
place with the precession frequency. Mass motions falling in these two categories 
should be avoided. 

6. An investigation of possible control techniques revealed that pure rate control 
and rate plus rate integral control would provide adequate damping of the spacecraft 
errors induced by internal disturbances. Initial attitude errors and attitude errors 
induced by external disturbances cannot be eliminated by these control techniques and 
require rate plus attitude control. 

7. Single-axis control was found to be acceptable for all control techniques and 
allows major reductions in the control gains and control system weight for near- 
cylindrical configurations spinning about a maximum axis of inertia. 

8. Optimization of the spin-up and extension technique for cable- or strut- 
connected spacecraft modules can lead to appreciable fuel savings for the extension and 
retraction process. Comparison of continuous -thrust, constant- rate, and constant- 
momentum extensions for an example spacecraft indicated that 22 lbm or 5 percent of 
the ideal spin-up and despin fuel could be saved by use of a continuous-thrust extension 
technique. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 17, 1967, 

125-19-01-20-23. 
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DEVELOPMENT OF THE LINEARIZED EQUATIONS OF MOTION 

The rotating spacecraft is considered to be the system of particles shown in fig- 
ure 1. A set of X,Y,Z axes fixed to the spacecraft is used to describe the rotational 
motion of the spacecraft with respect to a set of Xj, Yj, Zj axes which translate without 
rotation in inertial space and which remain parallel to a set of Xp.,Yjp,Zp axes fixed 
in inertial space. The general moment equation (ref. 6) about the origin of the X,Y,Z 
coordinate system is then 

5= Z f s*ai( m A) <A1) 

It is assumed that the system mass does not change during the time periods of interest 
so that 


M=^r j Xm j R j 

The absolute vector acceleration Rj is given by 

Rj = Rq + rj -}■ £2 x rj + 2 £2 x rj + £2 X ^£2 x rj^ 


(A2) 


(A3) 


and substitution of equation (A3) into equation (A2) yields 

M = pj x m jR 0 + ^ x m i(^ x r i) + 1 x m i (? x ? i) + Z r 3 x m 3 |p x ( s x r jj 


^r j Xm j (£2X^) + ^r j Xm j ?j 


(A4) 


or 


M=m s F s XR 0 + ^F j Xm i (nxr j ) + jF i xm i (Sxr,) + ^rjXmj(s5xf j ) 


X 5x [rj x mj(5x rjj] x *]) + £ 


r i x m i r i 


(A5) 


The acceleration of the origin Rq is found from the general force equation 
(ref. 6) __ 

P ~ ni s Rq + r s + £2 X r s + 2£2 X r s + £2 X ^£2 X rg^ 


(A6) 
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and the first term of equation (A5) may now be written as 

m s r s x Ro = r s x ^ “ r s x m s(« x r s) ~ r s x m s(^ x Tgj - r s X m s |f2 x (si x r s jj 

- r s X m s ^f2 s X r s ^ - r s X m s r s 

= r s X P - r s X m s (?2 X r s ) - r s x m s (f2 X r s ) - r s X m s p X f s J 

- S2 X j?g X m s (ft x r s jj - Q X m s ^r s xi^-r g X 


m s r s 


The equation of motion becomes 


M = r g X P + 1 


r i Xm 


- r s x m s . 


^ r j X m^Sl x rj^ - r s x m s ^2 X r s |j 
j(& X - r s X m s (& xr s | + r. X m^S2 X r^j 


+ < x 


^2 X r s J +S2X ^ ?j x mj^fi x rjj - r s x x r^j | 

[X m J ? J X " m s F s x *s] + [j 


r j X m i r j ■ r s x m s r s 


where the vector from the center of mass to the origin is 

mi 


- V __1 r • 
L m s 3 


(A7) 


(A8) 


(A9) 


To reduce equation (A8) to a more useful form, the particle system is represented 
as a large mass associated with the spacecraft and fixed with respect to the X,Y,Z axes 
and n smaller masses which move relative to the X,Y,Z axes. The rigid-body 
angular momentum vectors of the spacecraft, the n moving masses, and the spacecraft 
mass center are designated Hq, H n , and H s , respectively. The rigid-body angular 
momentum of the system H referred to the system center of mass is then 

H = H 0 + H n - H s (A10) 
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where 


Hg - r g X m s ^fi s X rgj 


n 

5 n - I F i x ”#j * F i) 

i=i 






Ho= 1 ?j * »#j X F i) 

j=n+l 


(All) 


and equation (A8) can be rewritten as 


dt 


r 

“ n 


"n _ 

1 

)q X 

Y m j ? j X ' m s*s x ' r s 

+ 

Y x m j*j - ? s x m s'r s 

1 

/ 


r 1 


1 


y. 




J 


(A 12) 


with 


n 

= Y^r. 

/_j m s 1 


j=l 

In component form, the equations of motion are as follows: 

M x = Ys^z “ z s^y + *x^x " *xy^y “ ^xz^z + *x^x ” ^xy^y “ ^xz^z 


- S2 z ^Iyf2y - Iy Z f2 Z - Iyx^xJ + £2y^. 


■z ai, z _ A zx' J “x ~ A zy^y) 


" Il7Y^*^Y “ I n 


Y m 3( X ^j ■ y 3 X j) 

i=i 


•Jit 


| Z m i( x f i - z f j) 


^ m l( y i 2 i " z 3 y j) 


m s( x s y s - y s k s) 


m s( x s z s ” z s x s) 


m s(y s ' z s - z s y s 


>% 


t Qny 


(A13) 


(A14) 
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• • * • • • 

ZgP x - x g P z + Iy^y - lyz^z " lyx^x + *y^y - %z^z “ lyx^x 

- - Izx^x “ J-zy^y) + ^z(lx^x ~ Ixy^y “ Ixz^z) ( 


'(’ , ‘‘"" 1 ) 


m. 


;(y S Xs - x s y s )| 


^ X 


J - z jyj) - m s(y s z s - z s y s )| n 


+ 



m c 


z s x s 



M z - x s^y " Ys-^x + *z^z ~ *zx^x " *zy^y + *z^z - *zx^x " ^zy^y 
~ ^y^x^x "* *xy^y " ^xz °z) + Dx^IyDy “ Iyz^z“ lyx^xj 


^ m i( z i x i " x i z i 
i=i 


m s ^s x s - x sZs) 




^ m i( z jyj - yjzj) 
j=i 


X m i{ x ih ~ y i x i) 

3-1 


m, 


i(z s ys - y s z s) 


m s(x s y s - y s x s )| 


1 



(A 15) 


(A 16) 
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where 


and 


-A 


*x ~ 1x0 + 


h = Vo + 


h ~ l zO + 


I m i( y i 2 + z i 2 ) - ms ( y 

i=i 


2 2 
s +z s, 


Z m i( x i 2 + z j 2 ) - ms ( x s 2 + z s 2 ) 
li=i 

2 m i( x i 2 + yj 2 )- m s( x s 2 +y s 2 ) 

j=i 


IxZ “ 


lyz - 


J xy = 


X m i( x i z i) ' m s( x s z s) 

li =1 

"n 1 

I m i(yj z j) * m s(ys z s) 
i=i 

"n 

X m i( x i y i) ■ m s( x s y s ) 

3=1 




II 

n 

£ m j ( y jyj + z j z j) • m s(y s y s + z sz s ) 

j=l 

“\ 


“n 


(M 

II 

2^ m^xjfcj + ZjZjj - m s (x s x s + z s z s j 

3=1 

“n " 




i z = 2 

Y m j (x j x j + YjYj) - m s (x s ks + Y s y s ) 



[3=1 J 

\ 

Ixz = 

~ n 

^ m j ( x j z 3 + z j x j) - m s( x s z s + z s x s) 
>1 


II 

SI 

_ n 

£ m j(yj z 3 + z jyj) - m s(y s z s + z s y s ) 


• 

*xy = 

j=i 

” n 

y m j( x jyj + y j x j) - m s( x s y s + ys x s) 



i=i 

J 


(A 17) 


(A18) 
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The coordinates of the jth moving mass are Xj,yj,Zj, and the coordinates of the mass cen- 
ter are 


n 




‘ s= 

j=l 



n 



(A19) 


For the special case when only one mass m with coordinates x,y,z is moving 
with respect to the spacecraft, the equations of motion yield 


Mx m s (y^ z ” zP y) + X x^x " ^xy^y " *xz^z + ix^x ” ixy^y “ ixz^z 

- f2 z ^Iyf2y - Iyz^z - Iyx^x) + fy^z^z " *zx n x " ^y^) 

+ Q [{xy - yx)Oy + (xz - zx)S2 z + (yz - zy)j 


(A20) 


M y m s ( zp x “ xP z) + ly^y “ lyz^z " *yx^x + iy^y “ *yz^z “ iyx^x 
fi x(l z ^z - *zx^x " *zy^y) + ^z(^x^x “ *xy^y " *xz^z) 

+ Q [(yx - xy)J2jj + (yz - zy)O z + (zx - xz)] (A21) 

and 

M z = ^(xPy - yP x ) + Iz^z ~ Izx^x ~ Izy^y + Iz^z ” Izx^x _ Izy^y 

- Oy^I x S2x - Ixy^y ~ Ixz^z) + ^x^Iy^y - lyz^z - Iyx^xj 
+ Q [(zx - xz)J2x + (zy - yz)S2y + (xy - yx)J (A22) 
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with 


and 


where Q is given by 


l x = ho + Q(y 2 + z2 ) 

fy 88 Vo + q (* 2 + z2 ) 

iz = i z o + Q( x2 + y 2 ) ^ 

I X z = Q(xz) 

Iyz = Q(yz) 

Ixy = Q(xy) 

i x = 2 Q(yy + zz) 

• • « 

Iy = 2Q(xx + zz) 

I z = 2Q(xx + yy) 

Ixz = Q(xz + zx) 
iyz = Q(yz + zy) 
ixy = Q(xy + yx)^ 

n _ m(m s - m) 


(A23) 


(A24) 


(A25) 


The spacecraft equations of motion (A14) to (A16) can be solved for the body rates 
fix, Oy, and O z . The motion of the rotating spacecraft is then defined in terms of the 
modified Euler angles 9, and cp. These angles, as shown in figure 2, relate the 
moving body axes X,Y,Z to the intermediate reference axes X^Y^Zj. From figure 2, 
the time derivatives of the Euler angles are 

cp - ^x + ^y tan 6 sin cp + O z tan 9 cos cp 

6 = fiy cos cp - i 2 Z sin cp ) (A26) 

\p = cos cp sec e + Qj sin cp sec 9 

The Euler angles found from equations (A26) and the body rates found from equa- 
tions (A14) to (A16) completely define the rotational motion of the spinning spacecraft. 
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For a large number of practical applications, motions involving small oscillations 
of the spacecraft spin axis from an equilibrium reference position are of primary interest. 
If it is assumed that the spacecraft spins about its Z-axis, that the Z-axis is initially an 
axis of maximum or minimum inertia, and that the Zj-axis is selected as the inertial 
reference, then 

sin <p =* tan cp ~ cp cos (p ~ 

sin 9 ~ tan 9 ~ 9 cos 9 ~ 



and 




Sly 


« fi. 


(A28) 


for the small oscillation regime. Consistent with these assumptions, the variable inertia 
terms can be considered to be sufficiently small in comparison with the spacecraft 
moments of inertia so that they may be neglected when multiplied by the oscillatory body 
rates or any angular accelerations. 


The reduction of the nonlinear governing equations to linear approximations can 
best be accomplished by first converting equations (A14) to (A16) to nondimensional form. 
As was done in reference 4, a nondimensional time r, an inertia term gpq, and a non- 
dimensional mass juj may be introduced by 


T _ ^zo* 

c _ T pq 
pq ~i z o 






(A29) 




mo 


where p and q range over x, y, and z. The remaining nondimensional terms are 
then 


Mr 


L p = 


^zO^zo 

Im2 


Uj = Xj C 


v j = yj\ 


°P = ~ 2 

“zo 


Pn = 


P fi Z o 2 PzO m s 

i- = k 3 .pis' 
j n zoPzO 


w i = z j\ 


€ nn ~~ 


A pq 


IzO^zo 


x i i m s 


J a zo 2 Vzo 


(A30) 
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and the nondimensional equations of motion become 


Lx = v sPz ~ ™ S Py + - e 


xy w y “ e xz U) z ~ e x a ’x “ e xy w y ~ e xz w z 


- w z ^eyCt)y - 6y Z o; z - £y X w x ^ + Vy^z^z - e zx w x “ e zy w y j 


n 

1 

\ ' / • • \ 

/ • * \ 

Z W V j ‘ V 3 U 3 
j=l 

- (u s v s - v s u s l 


n 


X Mj( u jWj - Wjhj) 

j =1 


+ J X ^(vjWj-W^,). 

r 1 


- (u s w s - w s u s j>w 


- ( v s™s - w s v s ) 


Ly - w s p x - u s p z + e y Wy - fyz^z " e yx^x + ^y w y " ^yz w z " ^yx w x 
" w x( £ z w z " e zx a ’x “ e zy a, y) + “z^x^x “ e xy w y ” e xz w zj 


n 

‘■N 

Z M i( v iVYi) 
li =1 . 



n 


X “ w ^) 

i=i 


+ | f - u fi) 


- ( v s w s - w s v 8 )>«z 


- (ws^s - U sW S ) 
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and 


L* = 


• • 


z ~ u sPy " v s Px + e z w z “ e zx w x - e zy a, y + e z w z “ e zx w x - e zy 0J y 
'y( e x w x " e xy co y ~ e xz w z) + w x( e y w y " e yz co z“ e yx w x) 


- ov 


+ < I - u rj ) 

Ir 


( w s u s - u s w s ) L } 


where 


n 


X Hi(wjV r v jWj) 

i=i 

Z Mj(ujVj - VjUj) 

j=l 


- (w S V S - V S W S )j>0)y 

~N 

- (u s v s - v s u s ) 


e x “ e XO + 


£y — £yo 


2 ^( V i 2 + W i 2 )- ( v s 2+w s 2 ) 


i =1 

" n 


I f‘j( u i 2 +w r)- (“s 2 +w s 2 ) 

j =1 



^X1 


e z e zo * 

2 "i( u i +T i 

) - (u s + v s 2 ) 





e xz 


e yz “ 


e xy - 


n 


X fl i( u J w i) - (% w s) 
3 =1 

~n 

Z Hj(VjWj) - (v s w s ) 

r 1 - 

~n 

I ^(Vi) - ( u s T s) 
1=1 


(A3 3) 


(A34) 
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2 2 

In accordance with the small oscillation assumptions, take u> k , 6<p, co p , 6 , <p , 

n n n n 

L p , p p , and the variable inertia terms Y pjUj 2 , N MjVj 2 , MjWj 2 , / MjUjVj, 

j=l 3=1 3=1 3=1 

n n 

^ PjUjWj, ^ PjVjWj, and their derivatives to be of order A. Here A is restricted 
3=1 3=1 

to be sufficiently small so that terms of higher order than A may be neglected in the 
governing equations. 

Thus, it follows that 

Lx = ysPz ~ w sPyj + ~ v e xy^Yj ~ y^xz^z^, + ~ ^xy^Yj ~ y£xz^zy 

o[a] oM °0 oM oM o[a3 oM o[a] 


w z ^y^y, ” ^yz^z y _ v e yx^x) y + w y “ v 6 zx^°x y ~ zy^y^ 

o[a] o[a] oM o[a] oM o[a 3 ] 


±x 

£ Mj(vi - T,Uj) 

J-l 


( u S v S - v s u s )i 


oM 



~n 

I 


\ / • • \ 



L ^3 ( U j W 3 - W j U j) 

- Y s w s - w s u s J' 

1 

3=1 


L 


J 


o[a] 



■* 

"N 


n 


\ 

Z ^( v iVYi) 

- (v s w s - w s v s ) 


i=i 






(A3 5) 


o[A] 
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o[a] o[a2] o[a] o[a 2 ] o[a 2 ] o[a 2 ] ofxf^ 

- “x( 6 z ^3 - ^zx^x,- ^_zy^y)^ '^cy^y/ ^xz^z) 

0[a] o[a3] qCa 3 ] OH o[a2] 


and 


1 - UjVj) 

j = l 


- ( v s"s - K s v s f> 


o[a 2 ] 
" n 

2 Mj(vjW 3 - WjVj) 
3=1 


- (v s w s - w s v s )L z 


O0 

" n 

I ^(wjUj - UjWj) 

3=1 


- (w S U S - UgWgj N 


(A36) 


Q[A] 

^ = v u sPy - v s P X/ + ,6z^z, - ,^A ,-.f zA,t .^z. - 

o[a] o[a!] o[a] o[a 2 ] o[a2] 0 [a] oM o[a2] 

' w y ^x^x, ' vfxy^y,' ^xz^z), + “x ^y^jr ■ v «yz‘‘'z,- «yz“x) 

o[a 2 ] 0 [a 3 ] oCa 2 ] OCA 2 ] 0 [a 2 ] o[a 3 ] 


2 Mj( w jVj - vjWj) 

j=l 


( W sV S - v s w s )i 


o[a2] 


2 Mj(wjiij - Uj Wj) 

3=1 ' 


( W s^s “ u s^s)J> w x 


V 


o|a2] 


I - v 3 “i) 

3=1 

C _ _ 

~°S 


- ( U sV S - VgUg) 


(A37) 
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or, if terms of higher order are neglected, that 

Lx = Ex^x + {^z ~ e y) w y a, z " w z(^xz “ e yz w z) 


2 ' w i“i) 

i' 1 


(u s w s - w s u s )\u> s 


n 


J Uj(vjWj - WjVj) 

j“l 


- (v s w s - w s v s V 


Ly - £yd)y - (e z - e x )w x a) z - w z (^yz + e xz w z) 


n 


I Mj (vjWj - WjVj) 

3=1 


n 


I M i( w i^i ' u j*j) 

i=i 


- ( v sW s - w S *s)>"z 


~ (ws'^s “ u S w s ) 


and 


= e z w z + ^z w z + 


n 


I Mj(uj'v j .- v jUj) 
]=1 


(u s v s - v s u s ) 


For the same range of disturbances, the Euler angles are given by 

= El) x + COy 8cp + ^ 


o[a] o(a5 oCa 1 /?] 


— = Wy 

dT v_x 


*z£, 

o[a] o[a1/2] 


+ 

dr 


z , T 


o[l] o[a 3 / 2 ] 




(A3 8) 


(A39) 


(A40) 


(A41) 
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and terms of order higher than A can be neglected, so that 


faT ~ w x + w z 6 

||=a)y-c o z cp) 

di// 
dr 




(A42) 


If the spacecraft dynamics are well conditioned, equations (A38) to (A40) and (A42) 
should give reasonable results in the small angle and rate regime, for which 

-15° § (p s 15° 

-15° s 0 < 15° 

or 

-0.0685 g A ^ 0.0685 

Solutions for higher values of A lead to correspondingly less accurate solutions. 

Since the form of the physical and nondimensional differential equations is identical, 
only the more convenient physical equations of motion are used. These equations can 
then be written as 

\ - h 


+ ~ \M X + fi z (i xz ~ ^yz^z) 


n 


I m i( z fj- x f 

i=i 


3 ( z 


z ij - m s (z s x s - x s z s ) 




1 m j( z fi - yfi)i 

j=i 


- m s(z s y s - ys z s)i ( A43 ) 


% - < ^ ~ \ M y + ^z (%z + *xz^z) 


T m j( z jyj - ypj) - m s (z s y s - y s z s ) 

Li— i 




X m i ( x i z i " z i x i) 

3=1 


- m s (x s z s - z s x s U 


(A44) 
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and 


]z 

Iz 




n 

] 

II 

isri*- 

M z + 

L 

2 m i( y J x i - x i y i) 
j=l 

- i»s(y s V - x sys)j 


(A45) 


where the moments of inertia I x and Iy, in general, may be approximated by their 
initial values in equations (A43) and (A44), so that 


*x ~ *xO + 


I- - Vo + 


n 


X m i(ym 2 + z 2 


U-i 

“n 


Z m i( : 

3=1 


r jO T "]0 


X. 2 . z. 2' 
]0 + Z J0 


m s(y S o 2 + z so 2 ) 


m s(x so 2 + z so 2 ) 


(A46) 


1 z *zO + 


X m i( x j 2 + yj 2 ) 

3=1 


L 2 . „ 2 


- m s lx s + y s 


J 


since the retention of variable inertia coefficients does not appear to add appreciably to 
the accuracy of the solutions. The other pertinent inertia terms are 



n 


Ixz ~ 

Z m i( x i z i) 

3=1 

- m s( x s z s) 


~ n 

r — v 


Iyz = 

Z m l( y i z i) 

j=l 

- m s(y s z s) 


(A47) 


J 


and the required time derivatives of the inertia terms become 


Ixz = 


!yz 


1-7 = 2 


I m ]{ x h + z i k i) 

j=i 

~n 

i=i 

~n 

X m j( x 3 x 3 + y 3^) " m s( x s x s + y S y S ) 
3=1 


- m s( x s z s + z s x s) 

~ m s(ys z s + z sYs) J> 


(A48) 
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APPENDIX A 


The Euler angle relations are given by 

cp - 

6 + Slgty = \ 

• | 

\f/ = fl z j 

If results should be required in nondimensional form, the transformations (A29) and 
(A30) can be applied to the solutions. 


(A49) 
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APPENDIX B 


PARTICULAR SOLUTION FUNCTIONS FOR THE UNCONTROLLED SPACECRAFT 


The forcing function for the differential equation with the constant moments of 
inertia approximation can be expressed by 

f 

F = £ AjFj(t) (Bl) 


j=l 


and thus the particular solutions are 


*\ 


(B2) 


F = 2 AjF,(t) 

i=i 

F = 2 AjFj(t)}» 

j = l 

f 

F = ^ A l F i« 

j=l > 

— s: £= 

where the functions Fj(t), Fj(t), and Fj(t) are given by the inverse Laplace transforms 


1 ' S 2 + X 2 / 


A 


f m _/ i \ f - 1 f F i (s) F ) (s> ■ 

' ' U : x2 j* ^2 + x 2 - prr?i 
fm ( 1 v 1 f sF i (s> sF > (s) ' 

J< '(o 2 -! 2 ) 1 is 2 + X 2 " s 2 + a 2 ! 




(B3) 


Several particular solution functions have been evaluated and are presented in tables I 
to EH. Other functions may be determined from equations (B3) if needed. 

The unit step function and the unit impulse function, which occur in these tables, 
are defined as 


U(t) = 


f 0 for t < 0 
1 for t > 0 


(B4) 
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and 


0 for t =* 0 





with 

6(t) dt = 1 


(B5) 


where 


t = 0 is the initial time of application of the step or impulse disturbance. 
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TABLE I.- SOLUTION TERMS FOR UNIT STEP AND UNIT IMPULSE FUNCTIONS 










TABLE III.- SOLUTION TERMS FOR THE TRIGONOMETRIC FUNCTIONS U(t)sin pt AND U(t)cos pt 


Fj(t) 


Fj(t) 


U(t)sin pt 

W*|p|*M 


-^[sinpt-Esin At] 


U(t)sin pt 

M-N 


sin Xt - t cos Xt] 
2X [X J 


U(t)sin pt 

|p| = M 


-^[sinat-2 sin xt] 


U(t)cos pt 
|A|#|p| #|a| 


„ U ^ O [cos pt - cos Xt 
X 2 - p 2L 


U(t)cos pt 

M=M 



U(t)cos pt 

M-M 


.[cos ot - cos At] 
X^ - a 2 
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TABLE IV.- FORCING FUNCTION COEFFICIENTS FOR CHARACTERISTIC PERIODIC MOTIONS 


Type of motion 


Motion 

coordinates 


Forcing function coefficients 


En 


Ei 


E 2 


e 3 


e 4 


Circumferential motion 
in plane perpendicular 
to Z-axis 


x = r Q cos pt 
y = r 0 sin pt 
z = z 0 U(t) 


(p + x y )(a+ p) 2 I rz 


(P + Xx)(g+P) 2 l rz 


(° + P) I 


■rz 


(p + v) ] 


rz 


lx 


rz 


Oscillation parallel 
to X-axis in 
XZ -plane 


x = r Q sin pt 
y = 0 i< 

z = z D U(t) 


[[a 2 + p 2 )p + (2crp)X^]l rz ) , ("[(o 2 + p 2 )x y + (2crp)p]l rz 1 


(2ap)I 


rz 


lx 


Plr 


Oscillation parallel 
to Z-axis in 
XZ-plane 


x = r 0 

v = 0 i 

{° 2 - P 2 Vrz 

1 

"(a 2 - p 2 )XyI rz 

j u L 

V J 

! 

Ix 

z = z 0 sin pt 


1 



Oscillation parallel 
to Y-axis in 
YZ-plane 


x = 0 

y = r 0 sin pt 
z = z 0 U(t) 


[(a 2 + p 2 )p + (2ffp)X y ]l 
ix 


[(a 2 + p 2 )x x + (2ap)pjl rz |^ 

(2ap)I rz 


Pl r z 


L J 


!x 


Oscillation parallel 
to Z-axis in 
YZ-plane 


x = 0 

Y = r o 
z = z 0 sin pt 


(a 2 - p 2 )pl r 


(a 2 - p 2 )Vr 


_L 


The term I rz denotes the constant inertia product Qr 0 z 0 . 
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TABLE V.- SOLUTION FUNCTION COEFFICIENTS FOR CHARACTERISTIC PERIODIC MOTIONS 


Type of motion 


B 4 


Solution function coefficients 


b 5 


b 7 


Circumferential 

hz 

( ( g + P) " ~ P^jj(p + ^y)(^y + ^x^x) 

^rz | 

J[(ff + p) 2 . (X 2 - p 2 |(p + Xy)(XIy - X x I x ) 

motion in plane 

2 >lxly 

l X 2 - P 2 

2XI X Iy 

l X 2 - P 2 

perpendicular 
to Z-axis 

► PX x Ix 

1 

■PVx^ 

J 

\ 


Oscillation 
parallel to 
X-axis in 
XZ-plane 


-ijj^ 2 + p 2 )p + (2 p)Xx](AIx + Vy) 


-i|[(o- 2 + p 2 )p + (2ffp)Xxj^XI x - Ayly) 


Irz( g + P? 2 ~| 
_ 2(x2 


T^[( p + Xx )i x + ( P + X ^ ^7^E P + Xx ) lx “( p + X ^ 


!_2(x 2 - p2)l x i : 
f (2a p) (ply 


y 

*x« 


(a2 + p 2 )(pl x + Xy Iy) 


2(x 2 - P 2 )lxly 
- (2ap)^ply - X x I x 


[(o 2 + P 2 )(pl x ' W) 


Oscillation 
parallel to 
Z-axis in 
XZ-plane 


IrzP(° 2 - P 2 )(XIx + XyIy) 


I rzP (a 2 - P 2 )(XI X - XyIy) 

i 

r rz( g2 - P 2 )(PI X + fyly) 

i 

- P 2 )(Plx - XyIy) 

2x(x 2 - P 2 )l x ly 


2X^X 2 - p 2 )l x Iy 


2(X 2 - P Z )I X Iy 


1 

to 

CO 

T3 


Oscillation 
parallel to 
Y-axis in 
YZ- plane 


{?•* + P 2 )p + (2ap)Xyj ^Xly + X x I x ^ 
- (x 2 - p 2 )pXIyj 


2X^X 2 - p2^I x Iy 


jj^ 2 + p 2 )p + (2ap)Xy] ^Xly - X x I x ) 

“(x 2 ' p2 W[:^X" 


2X(X 2 - p 2 )l xI y| 


2(x 2 - p2)i x iy] 

(2ap)^pl x + Xylyj 


"(° 2 + p 2 )(ply + X x I x ) 


2(x 2 - <J 2 )l x I, 

- (2ap)^pl x - Xylyj 


j^a 2 + P 2 )(ply - A X I X ) 


Oscillation 
parallel to 
Z-axis in 
YZ-plane 


I rz p ^ g2 ~ * 2 \**y * XxIx ) 


2x(x 2 - P 2 )l X Iy 


I rz p ( a2 " X 2 )(xly - X x I x j 


2x(x 2 - p 2 )l x Iy 


z(u 2 - P 2 ) (ply 


+ ^X^X 


2(x 2 - P^lxly 


[ rz(°^ ~ P 2 )(ply " *x r x 


2 (x 2 - p 2)] 


L x A y 


The term I rz denotes the constant inertia product Qr 0 z 0 . 




TABLE VI.- ASSUMED CHARACTERISTICS FOR MANNED 
ORBITAL RESEARCH LABORATORY 

[Conversion factors for SI Units are given on page 10] 


Parameter 

Module values 

Total 

values 

Manned 

Counterweight 

IxO slug-ft 2 

103 000 

30 000 

133 000 

IyO> slug-ft 2 

90 500 

73 000 

7 393 412 

W slu e - ft2 

173 000 

73 000 

7 475 912 

m s , slugs 

1 220 

557 

1 777 

Q, slugs 

36 


36 

a, rad/sec 



.4 



TABLE VII.- DISTURBANCE EFFECTS ON MANNED 
ORBITAL RESEARCH LABORATORY 
[Conversion factors for SI Units are given on page loj 




Error limits 

Disturbance 

Disturbance 

characteristics 

Rate 

rad/sec 

Attitude, rad 



Exact 

Approximate 

Exact 

Approximate 

Residual errors 

cp Q = e Q = 0.01 rad 
flxo = flyo = 0.01 rad/sec 

0.0161 

0.0161 

0.181 

0.182 

Step torques 

T x = 384 ft-lbf 
T y = 3300 ft-lbf 

0.0234 

0.0234 

0.0608 

0.0609 

Step inertia 

x Q = 50 ft 

0.00387 

0.00485 

0.0963 

0.0121 

products 

y G = 4 ft 

z 0 = 4 ft 





X-axis radial 

x = 45 + 5 sin ^ ft 

0.00591 

0.00686 

0.0102 

0.0107 

mass oscillation 

y o = 0 

z Q = 4 ft 





i 

Y-axis transverse 

x Q = 50 ft 

0.0196 

0.0224 

0.0720 

0.0857 

mass oscillation 

y = 4 sin ^ ft 
J 2 





. J 

z 0 = 4 ft 



. .J 


Vertical mass 

x Q = 50 ft 

0.0134 

0.0154 

0.0291 

0.0329 

oscillation 

y 0 = 4 ft 






z = 4 sin ^ ft 

mJk 
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TABLE Vm.~ ASSUMED CHARACTERISTICS FOR 150- FOOT SPACE STATION 


[Conversion factors for SI Units are given on page lo] 


Parameter 

Total value 

I xQ , slug-ft 2 

10 500 000 

I y0 , slug-ft 2 

10 500 000 

I z0 , slug-ft 2 

15 000 000 

m s , slugs 

2 270 

Q, slugs 

36 

a, rad/sec 

.314 
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TABLE IX.- DISTURBANCE EFFECTS ON 150- FOOT SPACE STATION 
[Conversion factors for SI Units are given on page lOj 


Disturbance 

Disturbance 
characteristic s 

Rate 

Error 

rad/sec 

limits 

Attitude, rad 



Exact 

Approximate 

Exact 

Approximate 

Residual errors 

(p Q = 0 O = 0.1 rad 

r = Ityo = 0.01 rad/sec 

0.0141 

0.0141 

0.176 

0.176 

Step torques 

T x = T y = 7500 ft-lbf 

0.0150 

0.0150 

0.0477 

0.0478 

Step inertia 
products 

x 0 = 49.5 ft 
y 0 = 49.5 ft 
z 0 = -20 ft 

0.00616 

0.00628 

0.0169 

0.0173 

Circumferential 
mass motion 

x = 70 cos ^ ft 

y = 70 sin ft 

J 70 

z 0 = -20 ft 

0.0153 

0.0162 

0.0380 

0.0518 

Radial mass 
oscillation 

x = 46 + 3.5 sin ft 

16 

y = 46 + 3.5 sin ft 

lb 

Zq = -20 ft 

0.00650 

0.00734 

0.0196 

0.0222 

Vertical mass 
oscillation 

x 0 = 49.5 ft 
y 0 = 49.5 ft 

z = -15 -5 sin ft 

5 

0.00590 

0.00724 

0.0175 

0.0199 


NASA -Langley, 1967 
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"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space ; The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 
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